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ABSTRACT 


This study addresses the problem of finding an energy conversion method to take 
advantage of the high maximum cycle temperatures achieved with solar central receivers. 
Most current practice is to use steam-based heat engines with solar receivers. Solar 
central receiver applications offer maximum cycle temperatures in excess of the operating 
range of steam-based Rankine cycles. Gas-based Brayton cycles are considered to take 
advantage of the expected increased cycle performance of higher temperatures. 

Simple Brayton cycles have been eliminated from consideration since they do not 
offer efficiency improvements over the lower temperature Rankine cycles. Derivatives 
of Brayton engines do offer favorable performance compared to Rankine cycles. 
Computer modelling was done to examine the effect of maximum temperature on 
efficiency of two Brayton cycle derivatives and a simple Brayton cycle (denoted as GT). 
The modified Brayton cycles include a combination of regeneration and one stage each of 
intercooling and reheat (IGT) and a Brayton cycle with steam injection (STIG). The 
turbine inlet temperature, the steam-to-air injection mass ratio (for the STIG) and the 
compression pressure ratios were treated as parameters in the analysis. The compression 
ratios investigated are 4, 8, 12, 16, and 20. Turbine inlet temperatures ranged from 
1000 K (1340'F) to 2500 K (4040°F). For the STIG the steam/air mass ratio varied 
between 0 and 0.5. Both First Law and Second Law efficiencies were examined as 
turbine inlet temperature and steam/air mass ratios were varied for each compression 
ratio. 

Results show that efficiencies were highest for the IGT followed by the STIG and 
GT, respectively. Considerable improvements in specific work output were 
demonstrated by the STIG over both the IGT and GT systems. First Law and Second 








Law analyses of the effect of turbine inlet temperature show a gradual increase of 
efficiency with turbine inlet temperature with diminishing returns at higher temperatures. 
These results demonstrate reason for optimism about increased performance for higher 
temperature air-based systems compared to Rankine systems at presently used 
temperatures. 
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NOMENCLATURE 


A Area 


b Specific availability (kJ/kg) 

C Constant 

Cp Specific heat at constant pressure (kJ/kg K) 

d Diameter 

GT Simple Brayton cycle 

h Specific enthalpy (kJ/kg) 

h Equations (7), (8) only -convective heat transfer coefficient 

HEX Heat exchanger 

hf Specific fluid enthalpy at saturation 

hfg Latent heat of vaporization (hg - hf) 

HRSG Heat recovery steam generator 

IGT Brayton cycle with intercooler, regeneration, & reheat 

m Mass flow rate 

Mr Steam/air mass ratio 

MW Molecular weight (kg/kmol) 

P, p Pressure (kPa) 

Pr Compression pressure ratio 

q, qjj Specific heat transfer (kJ/Kg) 

QHEX Specific heat transferred in the HRSG 

q sa t Specific heat transfer to heat water to saturation from ambient 

^superheat Specific heat transfer to heat water to superheat from ambient 

















































STIG 


Universal gas constant 
Total heat transfer resistance 
Specific entropy (kJ/kg K) 

Specific fluid entropy at saturation 

Difference between specific gaseous and fluid entropy at saturation 
( s g * s f) 

Brayton cycle with steam injection 
Temperature (K) 

HRSG approach temperature difference 
HRSG pinch point temperature difference 
Velocity 

Specific work (kJ/kg) 

Specific Carnot equivalent work of a heat input (kJ/kg) 

Quality (mass of vapor at saturation/total mass) 

Mole fraction 


Second Law efficiency 
First Law efficiency 
Dynamic viscosity 
Density 


Subscripts 


Numbers State points (see Figure 1) 
ap Approach point 

c Compressor 

p Pump 


































ACKNOWLEDGEMENTS 

I wish to thank my committee chairman Dr. Robert F. Boehm for his advice, 
guidance, and patience throughout the course of my work on this project. Appreciation 
is extended to my other supervisory committee members, Dr. James K.Strozier, and 
Dr. L. King Isaacson. I studied under all three professors. They made my stay at the 
University of Utah a rewarding experience. 








































INTRODUCTION 


The current practice for power generation from solar central receivers is to use 
steam-based heat engines. Even the newer receiver concepts, using molten salt or liquid 
sodium, are envisioned to use a heat exchanger external to the receiver to generate steam 
for driving a Rankine cycle system to produce power. In spite of this apparent 
contentment with steam-based systems, it is well-known that there should be a 
thermodynamic basis for expecting increasing cycle performance if the cycle maximum 
temperature is increased. However, any substantial increase in cycle maximum 
temperature may well preclude the use of steam. Gas-based cycles are therefore 
considered for this application. 

Air cycles have been considered for solar applications over the years. The Stirling 
cycle and Brayton cycle are both undergoing some development for solar power 
generation systems based upon dish (distributed) receivers. For larger systems such as 
might be used for solar central receivers, the development of Brayton-derivative cycles 
seems to be one of the few alternatives available. Brayton applications have been 
generally eliminated from serious consideration because they appear to offer too much of 
an efficiency penalty compared to lower temperature systems. This impression has been 
drawn because of increased receiver losses as temperatures are increased as well as the 
seemingly low cycle performance of simple Brayton cycles. However, there are some 
Brayton engine derivatives that may demonstrate favorable performance compared to 
Rankine systems. This study will evaluate the performance of two of these derivatives 
and compare them to the performance of the simple Brayton cycle (denoted as GT). 

























These include a familiar combination of a Brayton cycle with regeneration and one stage 
each of intercooling and reheat (denoted as IGT) and a less investigated Brayton cycle 
with steam injection (STIG). This study determines realistic performance of these cycles 
operating on air. Special concern is given to the impact of maximum cycle temperature 
on overall performance. 

As explained earlier, current solar power generation systems use a heat exchanger 
to transfer the heat from the molten salt or liquid sodium fluid to a steam based Rankine 
cycle. Part of the problem of making full use of higher temperatures afforded by solar 
systems is transfer of the solar energy to the power generation system. This study will 
not address the problem of how the exchange of heat from the initial fluid medium of the 
solar control receiver is made to the air or steam/air mixture working fluid of the cycles 
evaluated. This problem is a significant one and is not addressed here. In this study this 
heat exchanger is termed as a heater in the text and Figure 1. It will be assumed that the 
solar central receiver exists to heat the working fluid to the turbine inlet temperatures 
investigated in this study. 

The configurations of these three systems are illustrated in Figure 1. The 
advantages the IGT system enjoys compared to the simple GT system include the 
intercooler to decrease compressor work and the reheat to increase turbine work. Since 
both of these components incur an efficiency penalty, the IGT system includes a 
regenerator to increase efficiency by preheating the compressor discharge with the 
turbine exhaust. The performance characteristics of this Brayton derivative cycle have 
been well-documented. To serve as an example, several gas turbine texts were reviewed 
to compare the improvements offered by the IGT cycle over the GT cycle. Boyce [1] 
reports an 81.4% increase in specific net work and a 53.5% increase in efficiency at a 
turbine inlet temperature of 1256 K (1800°F) and a pressure ratio of 10.5 At a higher 
temperature of 1478 K (2200'F) and the same pressure ratio, he reports an increase of 
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52.5% in specific net work and a 56.2% increase in efficiency. Bathie [2] states that the 
IGT cycle provides a 33.3% increase in specific net work and a 42.2% increase in 
efficiency over the GT cycle at a turbine inlet temperature of 1400 K (2060T). Wilson 
[3] in his recent text lists performance curves for a Brayton cycle with intercooling and 
regeneration (no reheat). This system offered an increase of 20.1% in specific net work 
and an increase of 23.1% in efficiency over the GT cycle at a turbine inlet temperature of 
1465 K (2177°F) and a pressure ratio of 20. At a higher turbine inlet temperature of 
2051 K (3232°F) and the same pressure ratio, he reports an increase of 30.3% in 
specific net work and 41.4% in efficiency over the GT cycle. 

The STIG system also increases cycle efficiency by use of the turbine exhaust to 
generate steam. This steam is then injected into the working fluid, in this case at the 
heater, to increase the power output due to increased working fluid mass flow and also 
by increasing the working fluid's specific heat The theoretical performance of different 
variations of the STIG has been evaluated by several authors [4-17], It has been known 
for some time that water/steam injection increases efficiency and power output [4], 
Renewed interest started during the recent energy crisis period of a decade ago. Some of 
the earliest work on water/steam injection was performed by the NACA with the 
objective of augmenting the performance of gas turbine aircraft engines [5]. Early 
studies included water injection for power increase only, which typically decreased cycle 
efficiency 2% for each 1% increase in turbine mass flow [6]. The literature referenced 
here studies various other configurations of steam injection (5, 7-11] to include steam 
injection after the combustor [7], water injection by an evaporator prior to preheating in a 
heat exchanger [8], and steam generated by turbine exhaust which is used in both a gas 
turbine cycle and a steam turbine cycle [9]. 

The mentioned advantages of the STIG over the GT include the increases in cycle 
efficiency and power output. Boyce et al. [ 10] describe a 2-3% increase in efficiency 




















and a 12% increase in power output for each 5% increase by weight injection of steam in 
their external combustion steam injected gas turbine. Stochl [4] finds increases in 
efficiency and power output as high as 30% and 50%, respectively, at a turbine inlet 
temperature of 1366 K (2000^). Even more dramatic is Fraize and Kinney's [5] 
findings of efficiency increases of 8 to 14 points while increases for power output vary 
as much as 70% to 130% for their coal-fired gas turbine power cycle with steam 
injection. This also demonstrates that STIG applications may be possible for a variety of 
fuels and not just "clean" fuels. STIG systems have the advantage of simplicity when 
compared to the IGT or conventional combined cycles of gas and steam turbines. 
Existing gas turbine cycles have been modified for steam injection to handle higher loads 
with no increase in rotating machinery [12]. The STIG’s simplicity also means smaller 
capital costs in comparison to more complex systems Steam injection is known to 
reduce the level of pollutants when fuel is combusted. Field tests have found that a 5% 
by weight injection of steam will reduce the amount of NO x emissions to acceptable 
levels [10]. The STIG cycle has the disadvantages of corrosion and water consumption 
which all steam-based systems must cope with. Currently, two U.S. commercial firms 
offer STIG cogeneration systems. They are International Power Technology (TPT), Palo 
Alto, California and Mechanical Technology, Inc. (MTI), Latham, New York [13]. IPT 
has two operational installations of their Cheng-Cycle cogeneration STIG systems at San 
Jose State University and a Sunkist Growers installation in Ontario, California [13,14], 
There are no utilities currently operating steam injected gas turbines, but Pacific Gas and 
Electric Co. is interested in seeing the technology developed [ 13]. Actual testing of a 
STIG system in an electrical power generation configuration, with steam injected into a 
Westinghouse 191-6 gas turbine (compression ratio at 6.5, steam/air mass ratio at .05). 
yielded a 20% increase of generated electricity [15]. After 3000 hours of operation, 
there was no ill effect on the turbine blades. 







This study will vary turbine inlet temperature and steam/air mass ratio (for the 
STIG) for constant values of compression pressure ratio to evaluate their effect on cycle 
efficiency and power output. The compression ratios (denoted as Pr) investigated are 4, 
8, 12, 16, 20. This range is common to the studies in the literature. Brown and Cohn 
[16] state that the highest pressure ratio currently foreseen for single-shaft industrial gas 
turbines is 16. Although a pressure ratio of 20 would provide the highest efficiency, a 
pressure ratio of 16 was selected as the superior economic candidate for the basic STIG 
cycle. The turbine inlet temperature range investigated is 1000-2500 K (1340-4040°F). 
The temperatures examined in the literature were commonly in the range of 1073-1700 
K (1471-2600°F) with Day and Kidd [6] examining the range of 1228-2450 K 
(1750-3950°F). The investigated range is admitted to be high. The high end of which 
may not be technically possible without exotic materials and elaborate turbine cooling 
schemes, but an examination of the effect of high temperatures on performance is the 
purpose of this study. The range of steam/mass air ratios (denoted as Mr) investigated 
is 0 to 0.5. It will be shown here that maximum First Law efficiency is achieved within 
this range. This is confirmed by the literature. 

The Second Law efficiency will also be investigated for the configurations shown 
in Figure 1 plus two GT cogeneration applications where all the steam produced by 
turbine exhaust goes to process applications and is not injected in the gas turbine cycle. 
Second Law efficiency will be defined here as the ratio of the net work output 
determined by the First Law analysis and the maximum reversible work determined with 
the thermodynamic states calculated in the First Law analysis. This definition draws on 
the work of Boehm [18], Moran [19], and Kotas [20], Each of the systems to be 
investigated here, the GT, IGT, and STIG cycles, do not have combustion processes. 
These systems have heat addition processes from an external heat source. This source 
can be the solar central receiver and its initial working fluid. This study utilizes 










definitions from Boehm [18] and Moran [19] to express the heat addition processes over 
varying temperatures as reversible Carnot equivalent work. In addition it is necessary to 
complete the expression for reversible work by adding the availabilities of the inlet 
streams and subtracting the availabilities of the exit streams. The work of all three of 
these authors was referred to in expressing these availabilities. 






























CHAPTER n 


THEORY 

Computer modelling of a First Law analysis is the basis of this study. The model 
determines the thermodynamic states at each point in the three configurations. (The 
FORTRAN program GT developed for this study is in Appendix A). Given all the 
necessary states, the program calculates the performance characteristics of each 
configuration: specific work output, First Law efficiency, and Second Law efficiency. 
This chapter will detail the modelling techiques used in the computer program to 
determine the end states of various processes and phenomena. The assumed values used 
in the program are listed in Table 1. 

Compressor/T urbine 

The compression or expansion ratios and the assumed values of the component 
isentropic efficiencies are the basis of the calculation of the process (temperature change) 
through these components. They are assumed adiabatic. The assumed isentropic 
efficiencies are on the high side of those listed in the literature. For the IGT, the two 
stages of compression and expansion have the same pressure ratio which was 
determined by the square root of the total pressure ratio. The final turbine expands the 
working fluid to atmospheric pressure. The calculation for end temperature is shown 
below. 

Compm ssor 

T 2s = t i (Pr) k ' 1/k 


( 1 ) 











Parameter 


Assumed Value 


Ambient temperature, T j 

300 K 

Ambient pressure, 

101.325 kPa 

Air specific heat (initial iteration value), c p aQ 

1.0035 kJ/kg K [21] 

Air ratio of specific heats (initial iteration value), ^ 

1.397 [9] 

Air molecular weight, MW a 

28.97 kg/kmol [21] 

Air ideal gas constant, R a 

.28700 kJ/kg K [21] 

Steam specific heat (initial iteration value), c p s ^ 

1.8723 kJ/kg k [21] 

Steam ratio of specific heats (initial iteration value), 

1.327 [21] 

Steam molecular weight, MW § 

18.015 kg/kmol [21] 

Steam ideal gas constant, R s 

.46152 kJ/kg K [21] 

Water density, PHjO 

977.0 kg/m 3 [24] 

Universal gas constant, R 

8.31434 kJ/kmol K [21] 

Compressor efficiency, t] c 

0.87 

Turbine efficiency, Tl t 

0.89 

Pump efficiency, rip 

0.70 

Compressor pressure ratios, Pr 

4, 8, 12, 16, 20 

Pressure lost constant, kpj ost 

12.128 (kPa) 2 /K 
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Ti + T 2 ,-T. 

1 2actual ~ - 

Me 

(2) 

Turbine 

T 2 s = T i/(Pr) k ‘ 1/k 

(3) 

factual = ^1 ' ^t (^1 ■ ^2s)- 

(4) 


The specific work required and produced by the compressor and turbine 
processes, respectively, is calculated as 

w = c p factual' ^l) • (5) 

The program does the calculations in SI units. The specific work is in terms of 
kilojoules per kilogram of compressor air flow. The calculation of the turbine specific 
work as shown above is additionally multiplied by (0.95 + Mr). The 0.95 designates the 
95% of the air flow not lost from the working fluid for the turbine cooling (more on this 
topic later). The steam/air mass ratio, Mr, is added to the working fluid for the ST1G 
configuration. A variable specific heat for air or a steam/air mixure working fluid is 
determined using a variable air specific heat relationship [22] and a variable steam 
specific heat relationship [23]. These relationships use an average temperature of the 
start and end points. For the STIG cycle the mixture specific heat is determined using 
the steam/air mass ratio, Mr. 


Heat and Reheat 

The turbine inlet temperature is an independent parameter in this study so the end 
temperature of the heating process is known. The specific heat input required is 
calculated in similar fashion as the work calculation using an energy balance as shown 












































below. 


% ~ + Mr) Cp (T tur {jj ne j n i et - Tj) 


( 6 ) 


Heat Exchang ers 

The heat exchangers to be modelled in this study are the intercooler, the 
regenerator, and the heat recovery steam generator (HRSG). A simple but appropriate 
way was sought to express their efficiencies without setting any values for the heat 
exchanger's particular dimensions. A similar, analogous relationship was determined 
and applied to each heat exchanger. This was accomplished by expressing heat 
exchanger effectiveness in terms of the ratio of the minimum temperature difference 
between the two streams and the total convective heat resistance of the two streams. The 
conductive heat resistance of the heat exchanger is neglected. This ratio for each heat 
exchanger is set equal to the ratio of one representative heat exchanger. This is shown 
below. 


AT, 


mm 


R, 


AT, 


nun 


each 
I HEX 


representative 

HEX 


(7) 


This can be further expressed as 


AT 


min 


= R, 


AT, 


HEX 


min 


each 

HEX 


representative 

HEX 

















12 


—L + _JL 

hiAi h 2 Aj 


= AT, 


0110 1 + 1 


h 3 A 2 h 4 A 2 


For an assumption of equal heat transfer areas, the equation becomes 


representative 


-L + _L 

hi h- 


-L + _1 

h-i h A 


The representative minimum temperature difference selected was 50°C for the 
regenerator. A value of 51.1 *C was determined from a regenerator efficiency of 0.81 
[ 16] for a turbine inlet temperature of 1700 K and compression pressure ratio of 16. The 
representative values of the applicable convective heat transfer coefficients (Table 2) 
were selected [24, 25], The minimum temperature difference between the streams of the 
other two heat exchangers was calculated and is shown in Table 2. 

The locations of these AT m j n points are shown in Figure 2. The location is 
determined by the relative values of mass flow and specific heat between the two 
streams. These parameters determine the slope of the two streams (shown as linear on a 
T-length (x) diagram) and will determine where the two streams are closest. The 
locations of the intercooler and regenerator are pointed out but the location of the "pinch 
point" in the HRSG can be at various locations due to varying turbine exhaust 
temperature and steam/air mass ratio. This will be discussed later in this chapter. 


A method to express pressure losses without using any apparatus parameters and 
which was applicable at any location was sought. For steady flow through the com¬ 
ponents, the fluid velocity and geometry of the apparatus is generally constant. There- 
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Medium 

Convective Heat Transfer Coefficient (W/m^K) 

Forced convection of air 

200 

Forced convection of water 

2000 

Convection of boiling water 

5000 

Heat Exchanger 

Minimum Temperature Difference (*C) 

Regenerator (Air-Air) 

50 

Intercooler (Air-Water) 

27.5 


Regenerator (Air-Air) 
Intercooler (Air-Water) 
HRSG (Air-Boiling Water) 


26.0 




























































































fore, the Reynold's number is constant. 


Re = 


pVd 

-— = constant . 


( 10 ) 


For a constant Reynold's number and dimensions, the friction factor will be constant. 
The relationship for friction of internal turbulent flow is modified using an expression 
for mass flow rate and constant area. The following equation for mass flow rate is used. 


m = pAV . ( 11 ) 

The equation for the friction factor of internal turbulent flow is expressed as 

friction = f(Re) = = constant. (12) 

pV 2 

Applying Equation (11), Equation (12) becomes 


App 


= Ci . 


nr 


The ideal gas relation is used. 


_P 

Ap RT 


= C, . 


nT 


Since m and R are constant, the equation can be expressed as 


(Ap)p 

T 


= C 2 . 


( 13 ) 


( 14 ) 


( 15 ) 


This yields an expression for pressure drop in terms of pressure, temperature, and a 





































constant. A representative value of a pressure drop of 5 psi at atmospheric conditions 
was used to determine the constant, C 2 - The pressure drop through the intercooler, 
regenerator, heaters, and HRSG can now be determined using entry pressure and an 
average temperature. 

The literature used a number of techniques to model pressure losses in one or 
more of the major heat transfer components. Fraize and Kinney [5] set the pressure drop 
across both the combustor and steam generator at 2% of the initial pressure. Bhutani et 
al. [9] applied a pressure drop of 4% across only the heater. Boyle [11] accounted for a 
4% pressure drop across the combustor, a 4% drop in the turbine exhaust steam in the 
HRSG, and a 12% drop in the steam flow in the HRSG. He assigned generally larger 
values to account for parasitic pipe losses also. 

The method selected for this study to account for pressure losses has the 
advantages of determining the losses by use of flow parameters alone. This method is 
applicable to any of the mentioned components. 

Heat Recovery Steam Generator (HRSG) 

The effect of varying steam/air mass ratio and turbine exhaust temperature on the 
performance of the HRSG is constrained by a minimum pinch point limit (ATpp) 
between the two streams. The amount of heat to be recovered effectively by the 
incoming water requires that the temperatures of the exhaust gas stream must be higher 
than the water/steam stream temperature by a minimum amount at all points in the heat 
exchanger. Another temperature difference to be defined is the approach temperature 
difference (AT a p). This is the temperature difference between the entering hot exhaust 
stream and the exiting cooler water/steam stream. The location of the approach 
temperature difference as well as the possible locations of the pinch point temperature 
difference are shown in Figure 3. The literature describes the efficiency of the HRSG by 
setting AT a p at a certain value and then checks to see if the ATpp minimum limit has not 












Turbine inlet temperature 
is not hot enought to boil 
water in HRSG. 



Turbine inlet temperature 
is hot enough to start water 
to boil. Steam exits HRSG 
in saturation. 



Turbine inlet temperature 
is high enough for steam 
to reach superheat. 


re 3. Possible heat recovery steam generator (HRSG) pinch point (pp) locations. 
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been violated at the point where boiling of the water begins (point 6 in Figure 3). This 
sets the value of the exiting steam to the temperature of the turbine exhaust minus AT a p 
Once the temperature difference at the beginning of boiling approaches ATpp, the 
temperature of the exhaust stream at the pinch point is set at the water's saturation 
temperature plus ATpp. The temperature of exiting steam is then allowed to vary. 

The range of values of ATpp in the literature is 16.7 - 27.8*C (30 - 50*F) [4, 5, 7, 
9,11,21, 26]. Stochl [4] considers the value of 27.8’C (50*F) as a reasonable value for 
both AT a p and ATpp. In his study he examined four values for each. Other authors 
have set AT a p at a different value than ATpp including 41.7‘C (75*F) [16] and 55.5'C 
(100*F) [9]. The value of 26* C developed in this study for ATpp is certainly in line with 
the literature. In terms of the literature this study places both AT a p and ATpp at 26‘C. 
Figure 3 shows that the steam exiting the HRSG can be superheated, saturated, or can 
still be liquid depending on the turbine exhaust temperature and the steam/air mass ratio. 
The modelling of the HRSG in this study can handle each possibility. The modelling 
done in the literature commonly describes the performance of the HRSG only when 
^ap is set for exiting superheated steam and when the ATpp value is applied at the first 
point of boiling for steam exiting as superheated steam and saturated steam. In Figure 3 
these three situations are designated as pinch points 5, 6, and 3, respectively. It will be 
shown in the results of this study that the steam/air mass ratio which achieves the 
maximum First Law efficiency is the ratio which causes the pinch point (minimum 
temperature difference between the two streams) to move from the approach temperature 
location to the first point of boiling. This point of maximum First Law efficiency is 
confirmed by the literature. To describe the STIG performance over the entire range of 
turbine inlet temperatures and steam/air mass ratios the modelling of the HRSG was 
developed to handle all seven possible pinch point locations. 

The known conditions within the HRSG before calculations are made include the 
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turbine exhaust temperature, the ratio of steam to air in the turbine exhaust, inlet water 
temperature, and the water saturation temperature. The water is pumped to compressor 
discharge pressure plus the calculated pressure loss through the HRSG and combustor. 
The saturation temperature for that pressure at each of the selected compression ratios 
and the other saturation properties of latent heat of vaporation (hfg) and Sfg (for the 
Second Law analysis) were interpolated from Keenan et al. [27]. Additionally, a liquid 
water specific heat value for each compression ratio was determined by selecting the 
corresponding value of Cp for an average temperature (between ambient and saturation 
conditions) from the curve in Figure 2 in Keenan et al. [27]. The two unknown 
temperatures are the exiting steam temperature and the exiting hot stream exhaust 
temperature (denoted as T stac ^). These unknown temperatures are found by an iterative 
process. First the pinch point temperature difference (ATpp) of 26°C is set at first two 
possible locations (points 1 and 2 in Figure 3) and the unknown temperatures are found 
by equating the heat given up by the hot stream to the heat added to the cold stream. 
This assumes no heat loss from the HRSG. In both cases the temperature difference at 
the other possible pinch point location is checked to see if it is more applicable. If the 
heat transferred is sufficient to heat the water, for that stream/air mass ratio, to saturation 
then the next two possible pinch point locations (points 3 and 4 in Figure 3) for exiting 
saturated steam are checked. If the heat transferred in both of those situations is enough 
to heat the steam to superheat then the last three possible pinch point locations for this 
case (points 5,6, and 7, Figure 3) are investigated to see which is applicable. In general 
the specific heat (kJ/kg air) given up by the hot stream, which is equated to the total heat 
transferred, is calculated as shown below. 

^HEX = (0.95 + Mr) • Cp^^^e • (Turbine exit * ^stack ) (^) 

The specific heat added to the cold stream is determined as shown below (in the case of 
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superheat steam exiting). 


q - Mr [c pWater • (T sat - Tj) + h fg + c Rteam • ^exiting - T sat )] . 

steam 


(17) 


The HRSG mcxielling subroutine also finds the required heat input to raise the exiting 
steam to the desired turbine inlet temperature. This quantity is found by subtracting the 
heat transferred in the HRSG from the total heat required to raise the water from ambient 
conditions to the turbine inlet temperature. This is shown below. 

q| Hteam* Mr ,C Pwater <T sa<' T l ) + h fg + c P s team' turbine' T sat» ’ <)HEX 

Met (18, 


As mentioned earlier the values of specific heat at constant pressure (Cp) for the air, 
steam, and steam/air mixture streams are determined at average temperatures from tem¬ 
perature dependent relations for air and steam specific heat (Cp (T)) [22, 23]. 

Pump 

The required work to pump the water to the desired compressor discharge pres¬ 
sure is calculated as 


Specific Work = 


Mr • (Pr-Pj + APhrSG) 
Pwater • ^p 


(19) 


This quantity is negligible for the range of parameters of this study. 

Turbine Cooling 

For higher temperature applications in gas turbines a certain amount of compressor 
air is bled off from the working stream to cool the turbine blades to keep them within 
their materials' operating temperature limits. There are several modelling techniques 
listed in the literature. Fraize and Kinney [5] state that the amount of air bled from the 














compressor is 1.5% of the total compressor flow for each 10CFF above a turbine inlet 
temperature of 1750°F (1228 K). Day and Kidd 16] set a schedule of 4% flow loss for a 
turbine inlet temperature of 1500°F (1089 K), a 8.5% loss for a temperature of I750T 
(1228 K), and a 16% loss for a temperature of 200CTF (1367 K). For this study the 
technique of Bhutani et al. [9] was used. This is a straight 5% air flow loss to turbine 
cooling. This was reflected in earlier listed equations. This is an admittedly simple 
calculation. The authors of this technique applied it to their study which had a 1700 K 
(2600°F) turbine inlet temperatue. They make a good point that in reality the initial 
quantity of bleed flow would be larger but certain portions of the air coolant are regained 
to the working fluid at all subsequent turbine blade stages to the stage at which that air 
stream is applied for cooling. In keeping with the desire not to be locked in specifying 
any particular component parameters or dimensions than is necessary this technique of 
determining turbine cooling bleed off was used. 


The calculation of the First Law efficiency is a straightforward calculation of the 
ratio of net work output to heat input. For the IGT cycle there are two stages of 
compression and expansion as well as two heating components. The STIG cycle has a 
work requirement for the water pump and a heat input requirement to heat the steam 
exiting the HRSG to the turbine inlet temperature. The First Law efficiency calculation 


2 2 

j^l w turbine j ' £ j w compressor j' w pump 


j?l 9airj + 9steam 
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Second Law Efficiency 

At present there are no universally accepted definitions for Second Law efficiency. 
This study will follow the concepts of Boehm [18], Moran [19], and Kotas [20]. 
Second Law efficiency will be defined here as the ratio of the net work and the reversible 
work. The reversible work is further defined as the sum of the inlet stream availabilities 
and Carnot equivalent work of each of heat inputs minus the sum of the exit stream 
availabilities. 


w 


rev - ;^ ( b + w„-Ih 


inlets 


heat e exits 
inputs 


( 21 ) 


The results to be shown later will show the reversible work with and without accounting 
for the exit stream availabilities as a negative term in the reversible work. It may be 
argued that the exit streams that are vented to the stack are not available and need not be 
assessed as a loss or penalty. Some further definitions need to be made. First, 
thermodynamic availability is defined as the amount of work that can be produced when 
a substance at a given state reverts to the "dead" state while producing work with a 
Carnot cycle [18]. Here the dead state is usually defined as atmo r pheric conditions and 
the hypothetical Carnot cycle is assumed to be rejecting heat to the ambient temperature. 
T j. For this study the availability of a stream is expressed as showm below. 


b = h-h 1 -Tj (s - sj) . (22) 

The kinetic, chemical, and potential components of the availability are considered negli¬ 
gible. For a mixture of n ideal gases, the availability becomes (in molal basis) [19, 20] 



r T T i 

b= k=l y k [ 1 (1 ‘T )C Pk 

T 1 


dT] + R T j In (—) 


( 2 ?) 








The Carnot equivalent work of the heat inputs is calculated as 


T 1 

w = (l- T )q h • 


( 24 ) 


In comparison to this study's definition of Second Law efficiency, the definitions pro¬ 
vided in the literature are provided here. Moran [19] provides an example of a heat 
engine where work is produced by drawing heat from a high temperature reservoir, at 
temperature Tpj, and rejecting heat to a low temperature reservoir, at temperature T^. 
The Second Law efficiency can be expressed as follows (T Q is ambient temperature): 


W 


e = 


net 


(1 


T °)5Q h - (1 - “ )8 Ql 


(25) 


l H 


He states that this expression does not charge as a loss the availability exiting in the heat 
interaction at T^. This study follows more closely this alternate definition, which as 
he describes charges the availability exiting in the heat interaction at T^ as a loss. This is 
expressed as 


W 


e = 


net 


(I-^SQh 


(26) 


l H 


This study used this approach and assesses availabilities entering and leaving the con¬ 
trol volume of the system. Kotas [20] provides a Second Law efficiency definition for a 
gas turbine plant. This plant bums fuel internal to the cycle. His Second Law efficiency 
is expressed below. 
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e __ Rated plant net output _ 2 

Total availability of all input fuel components 

Kotas and El-Masri [28] provide details of the use of finding stream availability at each 
of the cycle state points to locate the cause and magnitude of losses within the system. 
This study uses a control volume approach to evaluate overall cycle performance, so 
finding the location and magnitude of the system's irreversibilities is not a goal of this 
study. This approach of a Second Law analysis is still quite new, but has a promising 
future. 

Simple Bravton Cycle (GT) 

For the GT and IGT cycles the referenced dead state is air at ambient temperature 
and pressure. The availability of the air inlet stream is then zero. The Second Law 
efficiency (e) for the GT cycle is defined as shown. 


w net 

w e a i r ~ b exhaust 


( 28 ) 


The Carnot equivalent work of the heat input is calculated as shown below. 


(• t 3* T! [ T 3 T, 

w Cajr =0.95j (1 -—) Sq-0.95 J (1 - -)c pair dT (29) 

T 

W 9aii = 0 95 >ir I(T 3 ' T 2>' T 1 ln tjr » ■ < 30 > 

The exhaust availability is calculated as 

b exhaust = °- 95 [ Ah - T l( As )]- < 3] ) 

♦Refer to Figure 1 to find the locations of referenced state points. 
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Since P 4 = Pj, then this term can be expressed as 


3 exhaust ~ ^-95 c p^j. ' ^i) - Tj in (^ )] 


As mentioned above, the availability of the inlet air stream is zero. In the IGT 
cycle there are two heating processes and the cooling of the air stream in the intercooler. 
These are expressed as three Carnot equivalent heat inputs. Each is calculated as shown 
above for the GT cycle. The exhaust gas availability is determined using the exhaust 
temperature from the regenerator. The cycle Second Law efficiency is defined as 


i?l We heatingj We intercooling Exhaust 


A Second Law efficiency for the STIG cycle is computed in an analogous manner 
to the approach covered for the GT and IGT cycles. In addition, two other Second Law 
efficiencies will be developed for comparison in which all the steam produced in the 
HRSG goes to process work elsewhere and is not injected into the gas turbine. In these 
second two cases, one will heat the exiting steam to the turbine inlet temperature and the 
other will not heat the steam beyond its HRSG exit temperature. 

In evaluating the Second Law efficiencies for these three systems, it will allow a 
comparison of a steam injectt i system and a cogeneration system. These systems are 
alike except the STIG system injects the steam produced. An evaluation of the level of 


irreversibilities of the steam injection can be made. 

For the first definition of the Second Law efficiency, the reference state is the 
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steam/air mixture at ambient temperature and pressure. In this case, the inlet air and 
water streams are at ambient pressure which is greater than their corresponding partial 
pressure in the exhaust mixture stream. The inlet stream availabilities are then nonzero. 
Using the earlier definition of availability (Equation (23)), the availability of the air inlet 
stream (on a mass basis) can be expressed as 

Tj * R In (1/y air) 

a ^ r inlet MWair 

The water inlet stream is expressed in a similar fashion. The availability of the mixture 
exhaust is expressed as 

Tctsrk 

^xhaus. = (Mr + °' 95) • C ftnix ‘ [T «** - T, - T, • In )]. (35) 

The heat input requirements for the STIG cycle are to the air and the steam. Depending 
on the turbine exhaust temperature and the steam/air mass ratio, the steam can exit the 
HRSG at superheat or in saturation. The Carnot work equivalent for the air heat input is 












The equivalent work for the heat input to take the steam through superheat to the desired 
turbine inlet temperature is calculated as 


"‘superheat = tb '' C Psteam [ ^ ' T *>' ) '' T > ' n ( ^ t » ' 
steam 


(39) 


The Second Law efficiency for this application is defined as 
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'sat steam 
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exhaust 


(40) 


In the last two Second Law efficiencies to be described for the GT cycle with stream 
production, the air and water streams never mix. The reference states are separate air 
and water streams at ambient conditions. The availability of these inlet streams is zero. 
The exit stream of steam is produced to do process work. Its availability will be added 
to w net in the numerator of the efficiency ratio. For this cycle, it was found that the 
steam could exit as liquid water, saturated steam, or superheated steam. This is due to 
the loss of increased turbine exhaust specific heat when steam is not injected into the 
working fluid. With only air in the hot stream of the HRSG, the HRSG modelling 
subroutine must be able to handle all 7 possible pinch points as described earlier. The 
availability of the HRSG exhaust stream is done as shown before except it is for only 
air. If the steam exits as liquid water its availability is described as 


1 water exit 

^water = c R va t er ^water * Tj) - Tj In ( ~ )] 

exit 


1 


(41) 


If the steam exists as saturated steam, the quality of the steam is found to properly ac- 




























count for the (h - hf) and (s - Sf) contributions to the availability by expressing them as 
x • hfg and x - Sfg. The quality is found as 

h-h f (h - hj) - (h f - h 1 ) qnEX * ^sat 


The availability is then expressed as 


^sat steam ^ { c pwater t^sat' ^i) * "^1 ( y )1 + x 1 hfg * Tj • x • Sfg}. 


For steam that exits as superheated steam (and similarly the case where the steam is 
heated to the turbine inlet temperature) the availability becomes the following relation. 


1 Sclt 

^superheat = ^ c Pwater l^sat * Tj) * Tj In ( ~y )] + hfg - Tj Sfg 
steam 1 


+ c Psteam ‘ l^exit * ^sat^ ' In ( M • (^4) 

*sat 

Additionally for the case where the steam is heated to turbine inlet temperature, the cal¬ 
culation of equivalent work of heat input must also handle the three possibilities of exit¬ 
ing steam (water, saturated steam, and superheated steam). The calculation for the case 
when the steam exits as liquid water will show how each case is done. 


T sat T 1 t 

^steam =*< + " ‘ ^ 
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The Second Law efficiency for the case of no steam injection and where the steam is 
heated to turbine inlet temperature is expressed as 


w net + ^steam(turbine inlet temp) 
We air + We steam' ^^exhaust 


(46) 


The Second Law efficiency for the case of no steam injection and where the steam is not 
heated beyond HRSG exit temperature is expressed as 

e _ w net + ^steam(HRSG exit temp) 

We air' ^ r exhaust 

This method of expressing the Second Law efficiency for these GT cycles with steam 
production closely follows the Second Law efficiency definition provided by El-Masri 
[28]. His definition is 


e = 
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net 


+ X b 
exits 


Qadded 


(48) 


Additionally he provides an interesting expression of the Second Law efficiency in terms 
of the Carnot efficiency and irreversibilities at each stage of the cycle. 


e ^Carnot • ^ L (49) 

components 


He expresses L as the ratio of the stream availability loss in a process to the total heat 
added to the cycle. 

This study presents several variations of the definition of Second Law efficiency. 
As a summary, these definitions are recapped in Table 3. Generally, each definition is a 
ratio of net work to reversible work. Depending on how the reference state is defined, 
the reversible work may include inlet stream availabilities. It will be shown in the next 




















































Second Law Efficiency Definitions 


Cycle Second Law Efficiency (e) a 


GT 


WorkOut ^ _ NetWork _ ncl 

Reversible Work " Carnot Equivalent Work of Heat Input w e . 


1GT 


_ Net Work _ 

Carnot Equivalent Work of Heat Inputs & Outputs 


w net 


X e heating ‘ c intercooling 
i=l 


STIG 


_ NetWork _ 

Availability of Inlet Streams and 
Equivalent Work of Heat Inputted 
to Steam and Air 


w net 

b air 4 b water + w e j ^ r 4 w c saturatw j 4 Wc S uperheat 
inlet inlet steam steam 


CT with Net Work and Sleam Availability _ w «l (rntrn aUel lemp) 

a, * We liquid * W Wa«d * "Wteattd 

Turbine water steam steam 

Inlet Temp) 


GT with 
HRSG 
(Steam at 
HRSG 
Exit Temp) 


Net Work and Steam Availability 

Equivalent Work of Heat 
Inputted to Air 


w net 4 b steam e»‘ tem P) 


"'air 


a Dcfiniiions do not include exhaust stream availability which will be shown to be not useful to this studs 
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chapter that including the exit stream availability in the reversible work is not useful to 
this study. The major components of the reversible work are the Carnot equivalent work 
of the heat added to the air and, depending on the cycle, added to the liquid water, 
saturated steam, and superheated steam. For the GT cycle with the HRSG to produce 
steam, the availability of the steam at either turbine inlet temperature or HRSG exit 
temperature is added to the net work. 


The limits to cycle performance are a result of maximum limits to the parameters of 
steam/air mass ratio and maximum cycle temperature. Both of these parameters are 
influenced by several variables. 


The appearance of condensation in the exhaust gas from the HRSG of the ST1G 
cycle places a maximum limit on the steam/air mass ratio. If the stack temperature at the 
HRSG exit is below the dew point temperature at ambient pressure then possibly 
corrosive condensation will occur within the HRSG. Although not the case in this 
study, when fuel containing sulfur is burned, any condensation in the exhaust flow will 
be acidic. The literature has selected a number of safe minimum stack temperature limits 
to preclude the formation of acidic condensation. Stochl [4] places a minimum T stac k 
limit in the range of 250-300°F (394-417 K) while Bhutani et al. [9] placed the minimum 
limit at 275°F (408 K) and Day and Kidd [6] selected 310°F (428 K) as the limit. 
Depending on the degree of water treatment of the feedwater for the STIG cycle, the 
condensation of the steam in the HRSG can have a corrosive potential for long term 
operation The effect on overall cycle performance to preclude any condensation within 
the STIG cycle will be shown later. 
















































Visible Plume 

Depending on the ambient conditions, a visible exhaust plume from condensing 
steam can form at a stack temperature above the prevailing dewpoint temperature. 
Although the plume is only steam, a visible plume is a sign of possible pollution. To 
preclude a visible plume, a lower limit than d( >ign conditions would be placed on the 
maximum steam/air mass ratio. In areas where large fluctuations of ambient conditions 
occur, it may be necessary to install sensing equipment to monitor ambient temperature 
and relative humidity so that the steam/air mass ratio can be adjusted to eliminate a visible 
plume. 

Compressor Surge Limit 

In fixed-geometry machines not specifically designed for steam injection the 
relative increase in mass flow through the turbine with steam injection can result in an 
increase in compressor back pressure and therefore an increase in the compressor 
pressure ratio. As the steam injection increases, the compressor pressure ratio increases 
toward its surge limit. The maximum increase in performance obtainable with steam 
injection then depends on the surge margin available to the particular compressor [4j. 
This study did not place any limits on steam/air mass ratio as no particular compressor 
surge limits were specified. 

Complete Oxygen Combustion 

This limit pertains to cases when fuel is burned in the combustor and is included 
here for information only. Complete oxygen combustion serves as a limit to both 
maximum cycle temperature and steam/air mass ratio. For a certain steam/air mass ratio 
there is no longer any excess oxygen to bum additional fuel to heat an increase in steam 
injection to the desired temperature. Further additions to steam are not possible if the 
turbine inlet temperature is not reduced [6,11]. 
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Materials' Temperature Limit 

As discussed earlier, turbine cooling is necessary for higher temperature 
applications in gas turbine cycles. Although over a billion dollars has been spent since 
World War II on metallurgical research in developing new alloys to extend jet engine and 
gas turbine temperature limits, there has only been an advance of approximately lOO'C 
[29]. The maximum permissible metal temperature in gas turbine blades has been 
extended from ~700 to ~ 80CTC (~1292 to ~1472‘F). Boyle [11] sets a maximum 
temperature limit at 1089 K (1500°F). This temperature limit is based on the fact that in 
40 years of research the same type of alloys of Fe - Cr - Ni are still primarily used. 
These alloys have a melting point around 1643 K (2550‘F) giving a sintering 
temperature (local welding at 75% melting temperature) of - 1230 K (~ 1750°F). Fraize 
and Kinney [5] state that at the time of their writing DOE programs were ongoing to 
extend ceramic technology to applications for gas turbines. This will be a difficult task. 
Ceramics have excellent resistance to thermal stresses but are very brittle. The addition 
of aluminum to ceramic compounds is being tried to improve ductility [1]. Successful 
ceramic compounds could possibly extend maximum temperatures to 1561 K (2350°F) 
[3], Application of coated carbon-carbon materials may extend temperatures to 1700 K 
(2600°F) [3], At this point in time, turbine inlet temperatures can be extended as far as 
the cooling schemes can keep the turbine blades within the permissible temperature 
limits. Fraas [29] states that with air cooling there is a point of diminishing returns 
where efficiency decreases due to increasing compressor air bleed off requirements to 
cool the turbine. Boyce [1] reports that this point of diminishing returns is 1533 K 
(2500*F). Fraas [29] does place hope in water cooling of gas turbine blades which may 
allow turbine inlet temperatures as high as 2273 K (2600"F). The water cooling will 
require pumping which in terms of cycle efficiency will be a much smaller cost than 
bleed off from the working fluid. It is therefore unlikely that the higher temperature 

















limits of this study (range is 1000-2500 K) will be applicable for some time yet. 


Steam Dissociation 

Another limit to the maximum cycle temperature is the degree of dissociation of 
steam at the increasing turbine inlet temperatures. Bhutani et al. [9] placed a limit of 
1 100'F (1000°F as a practical limit) on the superheat of steam in Rankine systems due to 
the decomposition of the steam. Fraas [29] describes the primary hazard is that steam 
increasingly dissociates at temperatures above ~ 315“C (600°F) which allows the 
hydrogen to diffuse through the metal walls and frees the oxygen to attack the metal 
walls. The literature referenced [4-17] has not placed any temperature or steam/air mass 
ratio limits for the STIG cycle due to this phenomenon. In some application of gas 
turbines dissociation can be valuable from the standpoint of increasing the turbine 
volume flow relative to the compressor volume flow [29]. 



















CHAPTER III 


RESULTS AND DISCUSSION 

This chapter will describe the results of the First Law and Second Law Analyses 
of the performance of the three Brayton cycle systems. 

First Law Analysis 

Figures 4 through 8 show plots of First Law cycle efficiency versus turbine inlet 
temperature for the three systems at each of the compressor ratios. Figures 9 through 13 
show the plots of specific net work versus turbine inlet temperature for each pressure 
ratio. A fourth curve (dash line) is shown on Figures 4 through 13 to depict the effect of 
limiting the steam/air mass ratio so that no condensation occurs in the STIG cycle. Each 
of the plots for the STIG cycle uses the steam/air mass ratios which achieve the 
maximum First Law efficiency. This includes the plots of specific net work which 
reflect the net work at the steam/air mass ratio which achieve the maximum efficiency. 
Each of the plots does not cover the entire turbine inlet temperature of 1000-2500 K. 
The plots were cut off just before the temperature where a steam/air mass ratio of 0.500 
no longer provided a maximum efficiency for the STIG cycle. Performance curves for 
the IGT cycle were not provided at the smallest pressure ratio, Pr = 4, since at lower 
turbine inlet temperatures, the 2 turbines cannot provide enough work to run the 2 
compressors. 

On comparing the 3 cycles in Figures 4 through 13 it can be seen that both the 
Brayton derivative cycles (IGT and STIG) provide improvements in both efficiency and 
net work output over the simple Brayton cycle (GT). The IGT cycle is superior than the 
other two cycles in efficiency over the entire temperature range. The STIG cycle 
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rst Law efficiency versus turbine inlet temperature 





















































rst Law efficiency versus turbine inlet temperature (Pr = 12), 











































'St Law efficiency versus turbine inlet temperature (Pr = 16). 


































































































versus turbine inlet temperature (Pr = 4). 




















































Figure 10. Specific net work versus turbine inlet temperature (Pr 








































































































'pecific net work versus turbine inlet temperature (Pr = 20). 
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provides a considerably greater improvement in net work output than the IGT cycle over 
much of the temperature range. As the pressure ratio increases the STIG cycle achieves 
a maximum efficiency at lower steam/air mass ratios. At the highest two pressure ratios, 
16 and 20, there appears to be a discontinuity in the specific net work plots for the STIG 
cycle. In this area of higher pressure ratio and lower turbine inlet temperature the 
maximum efficiency is achieved at a different HRSG pinch point location than the rest of 
the studied regime. This pinch point location occurs at higher steam/air mass ratios. The 
higher steam/air mass ratios yield higher net work outputs. This phenomenon will be 
covered later. 

The GT and IGT efficiencies continue to increase with temperature but do so with 
a much diminished return. The STIG cycle as mentioned achieves a maximum 
efficiency. The performance of the STIG cycle is dictated by the performance of its heat 
recovery steam generator (HRSG). For the temperature range studied here, as steam is 
introduced in the HRSG at small mass ratios the turbine exhaust stream is hot enough to 
heat the smaller quanities of steam relative to air, to a superheated vapor. The streams 
are closest at pinch point #5 as shown in Figure 3. The relative specific heat values (Cp) 
and mass flow ratios of the two flow streams are such that the slope of hot stream as 
seen in Figure 3 is quite shallow. Pinch point 5 is the appropriate location for the small 
mass ratios and the high temperatures covered in this study. As the mass flow rate 
increases the turbine exit temperature increases due to the increase in the turbine working 
fluid's specific heat. The exit temperature of the steam leaving the HRSG also increases 
since its temperature is fixed to the increasing turbine exit temperature by the minimum 
pinch point difference (ATpp). With the pinch point at this location, the efficiency 
increases with the mass ratio. In other words, the gain in the turbine work is greater 
than the penalty of heating more steam to the desired turbine inlet temperature. 

The phenomenon of increasing efficiency with steam/air mass ratio does not go on 




























indefinitely. The specific heat of the turbine exhaust increases with die mass ratio. The 
slope of the hot exhaust in the HRSG, as shown in Figure 3, decreases as the steam/air 
mass ratios increases. The difference in the temperature of the two streams finally 
reaches the ATpp value at the first point of boiling. This location was noted as pinch 
point location #6 in Figure 3. With further increases in the steam/air mass ratio the 
minimum temperature difference remains fixed at this location and the exiting steam 
temperature will then decrease as the mass ratio increases. From now on increases in the 
turbine work are less than the penalty of heating increasing amounts of steam to the 
desired turbine inlet temperature. After the pinch point location transition occurs the 
cycle efficiency decreases as mass ratio increases. Throughout much of the studied 
regime the mass ratio, which causes this transition to occur, provides the maximum 
efficiency plotted in Figures 4 through 8. This is verified by similiar plots in the 
literature [4, 6, 9, 19, 11, 16, 26], 

If the steam/air mass ratio is still increased at lower turbine inlet temperatures the 
exiting steam transitions from superheated vapor to saturated steam. This is denoted as 
pinch point location #3 in Figure 3. With further increases in mass ratio the pinch point 
location makes one last transition to location #4. Throughout the temperature range 
investigated the steam always exits as a saturated or superheated vapor. At higher 
turbine inlet temperatures the pinch point location will change from #6 to #7 before 
moving to position #4. 

At higher pressure ratios (16 and 20) and lower turbine inlet temperatures the 
maximum efficiency is not achieved at transition from pinch point location #5 to #6, but 
at the transition from #3 to #4. Efficiency increases in this regime when the pinch point 
is at the point of first boiling. In this area increases in turbine work with increasing mass 
ratio are greater than the penalties of heating the air and steam to the turbine inlet 
temperature. The penalty of heating the air was decreased by the higher compressor exit 






































temperatures at higher compressor pressure ratios. Figure 14 shows plots of efficiency 
versus steam/air mass ratio for increasing turbine inlet temperatures at constant pressure 
ratio (Pr = 16). It can be seen here that the maximum efficiency is reached at higher 
mass ratios (the transition of pinch point from location 3 to 4) at the lower temperatures. 
As the temperature increases, less and less gain is made in efficiency as mass ratio 
increases. The plots of efficiency versus mass ratio get very flat until the maximum 
efficiency is achieved at a much lower mass ratio (the transition from pinch point location 
#5 to #6). The result of this quick transition from a high mass ratio to a lower mass ratio 
is seen in Figures 12 and 13 where the steep drop in specific net work occurs. 

The effect of limiting the steam/air mass ratio to preclude any condensation in the 
STIG cycle is also shown in Figures 4 through 13. The steam/air mass ratio is allowed 
to increase until the HRSG stack temperature is just above the water dewpoint 
temperature for the steam's partial pressure. In general, the dewpoint temperature 
increases as the steam/air mass ratio increases or as the steam partial pressure increases. 
The performance curves of the STIG cycle are affected when a high steam/air mass ratio 
yields the maximum efficiency. This occurs at the higher temperatures for the lower 
pressure ratio curves and at the high and low temperature regions of the higher pressure 
ratios, 16 and 20. The STIG performance curves with the dewpoint temperature 
constraint were calculated by the FORTRAN program called CONDEN found in 
Appendix B. This program checks the STIG performance data files calculated by the 
main program GT. If the stack temperature at the maximum efficiency is below the 
dewpoint temperature, the steam/air mass ratio is decreased until the a stack temperature 
is achieved which is above the dewpoint temperature. Program CONDEN uses a 
relation for steam saturation temperature as a function of saturation pressure given by 
Irvine and Liley [23]. 

Selected performance data calculated by this study's program are compared to 























































data at the same parameters presented in Bhutani et al. [9] and Boyle [11]. Parameters of 
temperature inlet temperature, steam/air mass ratio, compression pressure ratio, 
isentropic component efficiencies, ambient conditions, and HRSG pinch point 
temperature difference, were matched to those of the literature and the performance char¬ 
acteristics were computed. Table 4 shows comparison of data with that of Bhutani et al. 
[9]. Table 5 shows comparison of data with that of Boyle [11], As can be seen there is 
close agreement between the study's efficiency data and that of the literature. There is 
more difference in the results of net specific work. These are due to different methods in 
pressure loss modeling and how the working fluid specific heat (Cp) is found. 

Second Law Analysis 

By conducting a Second Law analysis of the performance of the Brayton cycles 
further insight was sought into the performance characteristics of these cycles. Second 
Law efficiency expresses heat inputs as Carnot equivalent work requirements. This 
means that for a certain amount of heat, one could, at best, only utilize the Carnot 
equivalent of that heat input. This assesses smaller heat input penalties in calculating the 
Second Law efficiency. A Second Law analysis also accounts for the availability of inlet 
and exit streams. This can assess further penalty. In general. Second Law efficiencies 
are higher than First Law efficiencies. But the Second Law efficiencies do not mirror 
First Law efficiency in most of the cases studied here. 

Figure 15 shows Second Law efficiency versus turbine inlet temperature plots for 
the GT simple Brayton cycle at the three highest pressure ratios (12, 16, 20). Figure 16 
plots the relation of the same parameters for the GT cycle but this time the Second Law 
efficiency is calculated without the availability of the exit stream. As was mentioned in 
Chapter II, accounting for the availability of this stream may not have useful purpose as 
it is exhausted out the stack and is not really available. For the case of the GT cycle there 
is a noticeable difference in whether it is accounted for in the denominator or not. When 

































































Table 5 

Comparision of Selected Data with Boyle [11] 


Mr 

W b 
w net 


W c 
'"net 

rf 

0.0 

109.00 

.3338 

121 

.333 

0.157 

187.654 

.4472 

180 

.412 

0.170 

194.518 

.4442 

203 

.430 


fPressure ratio is 16. Turbine inlet temperature is 2000’F. 
‘’This study's data (W net - BTU/lbm air). 

“Boyle's data. 

“Percent difference from Boyle. 

















































































































































the exhaust stream availability is deducted from the penalty side of the Second Law 
efficiency, the efficiency continues to increase with turbine inlet temperature. When the 
exhaust stream availability is not accounted for there is a maximum Second Law 
efficiency within the studied temperature range. This is due to the effect of converting 
the heat input to a Carnot equivalent work term which asymptotically approaches a limit 
as temperature increases. This provides the insight that when expressing heat penalties 
as Carnot equivalent energy values, the GT cycle has a temperature limit for maximum 
efficiency. 

Figure 17 shows the same set of plots for the IGT cycle. Here the dashed lines 
are the Second Law efficiency plots calculated without accounting for the exhaust stream 
availability. There is not much difference between these plots, whether this availability 
is considered or not. This is due to the small quantity this availability represents. 
Because this system has a regenerator which takes advantage of the hot turbine exhaust 
to increase efficiency, the cycle's final exhaust temperature is much lower than the GT 
cycle at comparable conditions. The exhaust stream availability is then much lower for 
the IGT cycle. It is also seen that there is no maximum Second Law efficiency achieved 
for the IGT cycle over the studied temperature range. The effect of two heating stages 
does not provide the limiting effect of Carnot energy representation that was seen for one 
heating stage. For this cycle the Second Law efficiency continues to increase as 
temperature increases. 

The results of the STIG cycle are again more complex than the other two systems. 
It is also a matter of the performance of the heat recovery steam generator. Three 
definitions of the Second Law efficiency for steam producing cycles were calculated for 
comparison. The first efficiency represents the STIG cycle as depicted in Figure 1. The 
second two efficiencies represent steam produced in a HRSG without injection into the 
gas turbine. These two cases represent steam produced for cogeneration purposes. On 




















































































the first of these two cases, the steam is heated to the turbine inlet temperature. In the 
second case the steam is not heated after it leaves the HRSG. To handle these last two 
cases, the study’s program must calculate the turbine and HRSG processes with only air 
as the working fluid. 

To more easily represent what the Second Law analysis predicts about the STIG 
cycle performance, efficiency plots are calculated versus steam/air mass ratio at a number 
of turbine inlet temperatures. This shows the effect of mass ratio and turbine inlet 
temperature on performance. The effect of pressure ratio is not as significient as these 
other two parameters but will be discussed. First Figure 18 shows the effect of mass 
ratio on First Law efficiency for a range of temperatures at a pressure ratio of 16. It has 
been explained in detail earlier how the maximum efficiency is achieved. 

For the plots of Second Law efficiency for the STIG cycle it was determined to 
calculate the efficiency without accounting for the availability of the exit stream. Figure 
19 shows the difference between the two calculations. The dashed line is the efficiency 
without the exit stream availability accounted for. At low mass ratios it is such a large 
term compared to the heat penalties that it shows that it is better not to use any steam. 

Figure 20 shows the effect of mass ratio on the first definition of Second Law 
efficiency. This represents the Second Law efficiency for the STIG system depicted in 
Figure 1. It can be seen that a maximum efficiency is not achieved at lower 
temperatures. A maximum efficiency is later achieved as the turbine inlet temperature is 
increased. At lower turbine inlet temperatures the Second Law efficiency continues to 
increase when the pinch point is at locations 6, 3, or 4. The difference in the Second 
Law efficiency and the First Law efficiency in this regime is that for the Second Law 
analysis the heat penalty side is not increasing as fast as the net work side of the 
efficiency calculation. For the First Law calculation the opposite is true in this regime. 
The reason for the dropping off of the heat penalty for the Second Law efficiency can be 
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Figure 19. Second Law efficiency versus steam/air mass ratio for STIG cycle at 
turbine inlet temperature of 1500 K fPr = 16). 
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found in its method of calculation. The effect of increasing turbine inlet temperature is to 
decrease the heat penalty for heating the steam to superheat, when the steam exits the 
HRSG as saturated vapor, with pinch point locations at #3 and #4. Increasing turbine 
inlet temperature also decreases the heat penalty to heat the steam to the desired turbine 
inlet temperature, when the steam exits the HRSG as superheated vapor, with pinch 
point location at #6. This results in a wider temperature range than the First Law 
analysis, where the efficiency increased for pinch point locations at points 6 and 3. 
There is a limit to the effect of decreasing these heat penalties by increasing turbine inlet 
temperature. With increasing turbine inlet temperature, the transition to pinch point 
locations 6 and 3 happens at higher mass ratios. With a higher mass ratio, the effect 
previously mentioned is reversed by the higher quantity of steam. 

Increases in pressure ratio have a similar effect in increasing the temperature range 
where Second Law efficiency continues to increase as steam/air mass ratio. For 
increasing pressure ratio, the heat penalty to heat the exiting steam through superheat to 
the desired turbine inlet temperature decreases, plus the heating penalty of heating 
superheat steam decreases. In the case of pinch point #6, the exiting superheat steam 
increases in temperature as the pressure ratio increases. For pinch point locations #3 and 
#4, the saturation temperature at which the steam exits increases as the pressure ratio 
increases. These increases in exit temperature decrease the required steam heat inputs as 
mentioned above. 

In total, the analysis of the Second Law efficiency shows that at lower turbine inlet 
temperatures the efficiency increases throughout the steam/air mass ratio range. At 
progressively higher temperatures, as pressure ratio increases, a maximum Second Law 
efficiency is finally achieved at the same steam/air mass ratio as the maximum First Law 
efficiency. 

For comparison purposes, a Second Law analysis is done on a GT cycle which 



















produces steam for cogeneration purposes but has no steam injection. Figure 21 shows 
similar Second Law efficiency versus steam/air mass ratio plots at various turbine inlet 
temperatures for the system which heats the steam to the turbine inlet temperatures. This 
figure shows that the Second Law efficiency continues to increase throughout the mass 
latio range for all the selected temperature plots. A maximum efficiency is reached at a 
higher turbine inlet temperature (not shown in Figure 21) than any of the other systems 
evaluated. 

This system experiences the same decreases in some of the heat penalties at pinch 
point locations #6 and #3 as explained for the previous systems. In addition, this system 
accounts for the availability of the steam at the turbine inlet temperature on the net work 
side of the efficiency ratio. This has a positive effect on the Second Law efficiency as 
turbine inlet temperature increases. By accounting for the steam in this manner, the 
Second Law efficiency continues to increase with mass ratio over a wider temperature 
range than the STIG cycle and the system to be shown next. 

Figure 22 shows the Second Law efficiency versus steam/air mass ratio plots for 
various temperatures at Pr = 16 for the system in which the steam is not heated after 
exiting the HRSG. For this system there is always a maximum Second Law efficiency 
achieved within the mass ratio and temperature ranges of this study. At lower turbine 
inlet temperatures, the maximum efficiency occurs at the transition from pinch point #3 
to #4. At this point the availability of the steam exiting the HRSG decreases with the 
pinch point located at position #4. At higher turbine inlet temperatures the maximum 
efficiency occurs at lower mass ratios where the pinch point location transitions from #5 
to #6. The system previous to this showed that higher efficiencies can be achieved by 
further heating the steam from the point left off at by this system. 

As a note, in calculating the performance of a HRSG with only air in the hot 
stream, the exiting steam can leave the HRSG as liquid water, saturated steam, or 












































temperatures (Pr = 16). 
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HRSG where steam is not heated beyond the HRSG at five turbine 
inlet temperatures (Pr - 16). 

















































































CHAPTER IV 


CONCLUSIONS AND RECOMMENDATIONS 

From the First Law analysis it can be concluded that both of the Brayton derivative 
cycles, IGT and STTG, offer significant improvements over the simple GT Brayton 
cycle. In this study a major emphasis has been placed on examining cycle performance 
as a function of cycle maximum temperature. No attempts are made to define equipment 
concepts for the heat addition process. As turbine inlet temperature is increased, all three 
cycles demonstrate a diminishing increase of efficiency as temperature increases but the 
IGT and STIG are quite superior to the GT cycle. The IGT cycle is the clear winner in 
First Law efficiency with the STIG in second. The STIG cycle shows a most dramatic 
increase in net work output. It demonstrates an increasing relation with temperature. 
The STIG work output can be more easily controlled by the steam injection rate, which 
need not require fluctuations in rotating machinery speed. The STIG cycle represents a 
much simpler system. This equates to smaller capital costs. The STIG cycle results 
shown here and the other advantages mentioned in Chapter I demonstrate that the STIG 
cycle is the most promising candidate for higher temperature applications. Table 6 
shows the improvement in the net work output and efficiency that the IGT and STIG 
cycles have over the GT cycle at selected maximum temperatures and a constant pressure 
ratio. 

As pointed out in Chapter II the limits on turbine inlet temperature due, to the 
effectiveness of turbine cooling schemes and the turbine material's temperature limit, do 
not allow application of the whole range of the temperatures investigated here. Further 
technological advances will have to be made before the application of the higher 
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temperatures studied here can be realized. Within current temperature limits the STIG 
cycle still exhibits a greater net work capacity although its greatest improvements over 
the IGT cycle in net work are achieved at higher temperatures. The STIG is not as 
efficient as the IGT in this range but offers several other improvements of simplicity, 
lower capital costs, and flexibility in the control of the steam injection to vary work 
output or turbine inlet temperature. These advantages are deemed more valuable overall 
than the IGT's improvement in efficiency. And with further advances in maximum 
turbine inlet temperature, the STIG cycle will offer greater improvements in net work. 

The Second Law analysis provided some further insights in the performance of 
these cycles. When the efficiency is represented as the ratio of net work output to 
reversible work the GT cycle experienced maximum efficiencies at intermediate turbine 
inlet temperatures. The Second Law efficiencies calculated for the IGT cycle was similar 
to its First Law efficiencies. Its Second Law efficiency continued to increase with 
turbine inlet temperature. Like the First Law efficiency, the increase in the Second Law 
efficiency diminished in the higher temperature range. The maximum Second Law 
efficiency for the STIG cycle was achieved at a higher steam/air mass ratio than the First 
Law efficiency at the lower temperature ranges for each pressure ratio. At the higher 
pressure ratios no maximum efficiency was achieved within the steam/air mass ratio 
range at these lower temperatures. This demonstrates that for the STIG cycle, less 
irreversibilities are experienced at higher mass ratios than those which achieved the 
maximum First Law efficiencies at these lower temperatures. Location of the higher 
irreversibilities within the cycle was not included in the scope of this study. 

Both First and Second Law efficiencies have been provided for each cycle and for 
a GT cycle with a HRSG. Comparisons are made of when maximum efficiencies are 
achieved as the parameters of turbine inlet temperature, steam/air mass ratio, and 
compressor pressure ratio are varied. But for solar power generation applications which 
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is more applicable, First or Second Law efficiency? The utility of the long used First 
Law efficiency is its value of measuring what you are getting for your money. The First 
Law efficiency tells the plant manager what amount of power is generated for the 
associated fuel costs. For a solar powered system there are no fuel costs. The Second 
Law efficiency makes more sense in telling the designer how close a given design is to a 
completely reversible system. And with a more detailed Second Law analysis of each 
system component the designer can locate where the major irreversibilities are. 

Where the parameters of this study provide a maximum in First Law efficiency but 
not in Second Law efficiency, it would be more useful to solar power applications to 
seek the parameters which provide the maximum Second Law efficiency as shown in 
this study. 

Each of the simplifying assumptions of this study represents recommended areas 
of further study in this area. There are two recommendations which should be among a 
more detailed study of these power cycles. First, a turbine cooling scheme should be 
used that varies the amount of compressor air bleed off with the turbine inlet 
temperature. Bhutani et al. [9] made a good case for the straight 5% bleed off of 
compressor air at a turbine inlet temperature of 1700 K. This study used this method. 
To provide a conservative performance analysis for higher temperatures, a cooling 
system which tasks more air for turbine cooling may be better. 

A second recommendation for further study is to adopt a set pressure loss figure 
for each applicable process from among the literature. Table 7 lists the pressure losses 
calculated for the STIG and IGT cycles at various parameters. The method developed in 
this study for calculating the pressure loss in terms of average temperature and initial 
pressure is applicable for only selective ranges of these two parameters. This judgement 
is made by comparing the results in Table 7 to selected pressure loss figures from the 
literature listed in Table 8. 
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Pressure Losses Calculated in This Study a 


Cycle 

Component 

Pr 

%Apb 

STIG C 

Heater 

12 

0.9 



16 

0.5 


HRSG 

12 

0.7 



16 

0.4 

IGT 

Intercooler 

12 

3.6 



16 

2.8 


Regenerator 

12 

0.7 



16 

0.4 


First Heating 

12 

1.2 



16 

0.6 


Second Heating 

12 

13.9 



16 

10.1 


^Turbine inlet temperature for each example is 1500 K. 

b%AP represents percent pressure drop (%AP = APUOOyPjnjtjai). 

c Mr for STIG cycle examples is 0.250. 









































Table 8 

Pressure Loss Schemes Listed in Literature 


Reference 

%AP a 

Comments 

Fraize and Kinney [5] 

2.0 

Assessed at each heating process and HEX. 

Bhutani et al. [9] 

4.0 

Assessed at combustor only. 

Boyle [ll] b 


Assessed at burner 



Assesed to turbine exhaust flow in HRSG 




Assessed to steam flow in HRSG 






























































There are numerous further adaptations and more detailed modeling schemes 
which are applicable to more advanced study in this area. Studying the effect of changes 
in ambient conditions would not require any adaptation to the existing program. A 
further change to the STIG cycle without an addition of rotating machinery is 
intercooling to reduce compressor work. Where this study left off in its Second Law 
analysis, further studies can analyze individual component irreversibilities. There is 
certainly more to be done in this area. A goal of this study is to present its findings in a 
manner to give further work in this area several possible stepping off points. 





































































PROGRAM GT 


Program finds 1st and 2d law cycle efficiency and net work versus 
turbine inlet temp for simple Brayton cycle (GT); Brayton cycle 
with intercooling , reheat, and regeneration (IGT); and Brayton 
cycle with steam injection (STIG). 

CHARACTERS CYCLE, OUTPUT 
INTEGER C,I,J,K,K1,L,LAW,II,JJ,R 

REAL T6,T(9),P(9),Te(9),Pe(9),T5s,T5se,ErA,EPS(6),Wc(2),Wct, 

+Wt(2),Wtt,Wp,Wpe,Wee,Wte,Pro(5),PL,Mr,Cpmix,Cph,Cpsc,Qha(2),Qhs, 
+Qh,Qhex,Qhat,Qhe,Qhse,Qhexe,Qhat e,Cpmixe,Cphe,Cpsce,Ke,Tstack, 
+Tstace,MAX(1334,5),MAX1(491,4),Wnet,Wnete,Be,We(3),Wt1,Bs1 ,Bs2, 

+Be1,Wes(3),T5s1,Tstak1,Qhexa,T9a 

DATA Pro(1)/4•0/,Ke/.4299/,SL41 /4•2857142E-04/,SL42/4.3103448E-04/, 
+SL81 /4-1121495&-04/, SL82/3.960396E-04/, SL121 /3- 8709677E-04/, 

+SL122/3.7719298E-04/, SL161/3-5251798E-04/, SL162/3.6885245E-04/, 
+SL201 /3.8E-04/,SL202/3*8636363E-04/ 

Make choice of which cycle to be analyzed. 

¥RITE(*,10) 

FGRMAT('0','Enter code below for desired cycle analysis:'/' ',T5, 

+'G-For simple Brayton cycle'/' ',T5,'I-Eor improved Brayton cycle' 
+,' with intercooling,regeneration,&reheat'/' ',T5,'S-Por Brayton ' 
+,'cycle with steam injection') 

READ(* f 20) CYCLE 
FORMAT(A1) 

Make choice of which compression pressure ratio to be analyzed. 
WRITE(*,30) 

FORMAT('O'.'Enter integer below for desired compressor pressure ', 
+'ratio (Pr):'/' ’,T5,'1-Pr=4'/' ’ T5,’2-Pr=8'/' •,T5,'3-Pr=12'/ 

+' *,T5**4—Pr=16'/' , ,T5,'5-Pr=20') 

REAL(*,40) L 
FQRMAT(II) 

Make choice of which type of output is desired. 

WRITE(*,50) 

FGRMAT('O','Enter responce for desired output below:'/' ',T5, 
+'T-Tabular data output'/' ',T5,'1-Plot of turbine inlet ', 
+'temperature vs efficiency'/' ',T5,'2-Plot of steam/air mass ', 

+'ratio vs efficiency'/' ',T5,'5-Plot of First Law efficiency ', 
+'and corresponding mass ratio'/' ',T7,'versus turbine inlet ', 

+'temperature') 

READ(*,60) OUTPUT 
FORMAT(A1) 

Make choice of which type of efficiency to be outputted. 













































IF(OUTPUT.EQ.’1 1 .OR.OUTPUT.EQ.'2') THEN 
WRITE(*,61) 

FCRMAT('O','For efficiency plots specify which type of 
+'analysis below:'/' ',T5,'1-First Law efficiency'/' ',T5,’2-', 
+'Second Law efficiency') 

READ(*,62) LAW 
F0RMAT(I1 ) 

MDIF 

Choose which turbine inlet temperature. 

IF(0UTPUT.NE.'1’.AND.OUTPUT.NE.'3') THEN 
WRITE(*,63) 

F0RMAT('O','Enter value of turbine inlet temperature(K)', 

+ '(F6.1).') 

KEAD(*,65) T6 
FORMAT(F6.1 ) 

ENDIF 

Choose which steam/air mass ratio (Mr), (for STIG cycle). 
IF(CYCLE.EQ.’S') THEN 

IF(OUTPUT.EQ.'T' .OR.OUTPUT.EQ. '1 ') THEN 
WRITE(*,70) 

FORMAT('O','Enter value of steam/air mass ratio (Mr)', 

+'(F5.3).’) 

READ(*,80) Mr 
FORMAT(F5-3) 

ENDIF 

ENDIF 

IF(CYCLE.NE.'S') Mr=0.0 
IF(OUTPUT.NE.'T') GO TO 230 

For tabular data output for specified turbine inlet temperature 
and specified Mr. 

CALL CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,Wt,Wp,Wnet,Cph,Cpsc, 
-+Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wt1,Bia,Biw,Be1 ,Wes,K1, 
+T5s1,Tstak1,Qhexa,T9a»Bs1 ,Bs2) 

DO 90 1=1,9 

Te(I)=1.8*T(I) 

Pe(I)=0.14504*P(I) 

CONTINUE 
Wct=Wc(1)+Wc(2) 

Wtt=Wt(1 )+Wt(2) 

Qhat=Qha(1)+Qha(2) 

Qh=Qha(1 )+Qha(2)+Qhs 

Wce=Ke*Wct 

Wpe=Ke*Wp 

Wte=Ke*Wtt 

Wnete=Ke*Wnet 

Qhate=Ke*Qhat 

Qhse=Ke*Qhs 





























































Qhe=Ke*Qh 
Qhexe=Ke*Qhex 
Cpsce=.2388*Cpsc 
Cpmixe=.2388*Cpmix 
Cphe=.2388*Cpb 
T5se=1.8*T5s 
Tstace=1 .8*Tstaek 

IF(CYCLE.NE.'I') THEN 
IF(CYCLE.EQ.'G') THEN 
WRITE(15,100) PL 

FORMAT ('O', 'SIMPLE BRAYTON CYCLE WITH OVERALL Pr= \F4.1) 
GO TO 120 
ENDIF 

WRITE(15» 110) PL,Mr 

FORMAT('O','BRAYTON CYCLE WITH STEAM INJECTION, Pr= \F4.1, 
+2X,'Mr= ',F5-3) 

WRITE(15,130) 

FORMAT('O','STAGE', 4X,'T(K)',6X,’T(R)’,5X,'P(kPa)',3X, 
+'P(psi)'/'O') 

1=1 

WRITE(15,140) I,T(1 ),Te(1 ),P(1 ),Pe(1) 

F0RMAT(' ',2X,I1,4X,F8.3,2X,F8.3,2X,F8.3,2X,F8.3) 

1=2 

WRITE(15,140) I,T(5),Te(5),P(5),Pe(5) 

IF(CYCLE.EQ.'S') THEN 
WRITE(15,150) T5s,T5se 

FORMAT(' ','STEAM IS INJECTED INTO COMBUSTOR AT \F6.1, 
+'K(',F6.1,'R)') 

ENDIF 

1=3 

WRITE(15,140) I,T(6),Te(6),P(6),Pe(6) 

1=4 

WRITE(15,140) I,T(9),Te(9),P(9),Pe(9) 

IF(CYCLE.EQ.'S') THEN 

WRITE(15,160) Tstack,Tstace,K,K1 ,T5s1 ,Tstak1 ,Qhexa,T9a 
FORMAT(' ','EXHAUST STACK TEMPERATURE IS ',F6.1,'K(', 
+F6.1,'R)'/' ','HRSGmode is ',11/' ','Air HRSG Results are:', 

+' Mode is ',11/' ',22X,'Water leaves at temp ',F6-1/' ',22X, 
+'Tstack is ',F6.1/’ \22X,'Qhex is ',F8.3/' ',22X,'Exhaust ', 

+'enters HRSG at ',F6.1) 

BIDIF 

WRITE(15,170) 

FORMAT('O'/’O',10X,'(kJ/kg)’,IX,'(BTU/lbm)') 

WRITE(15,180) Wet,Wee,Wp,Wpe,Wtt,Wte,Wnet,Wnete,Qhat,Qhate, 
+Qhs,Qhse,Qh,Qhe,Qhex,Qhexe,Be,We(1),We(2),We(3) 

FORMAT(' ','Wc',8X,F7.3,2X,F7.3/' ','Wp',8X,F7.3,2X,F7.3/' 
+'Wt',7X,F8.3,2X,F7.3/' ','Wnet',5X,F8.3,2X,F7.3/’ ','QHair', 
+4X,F8.3,2X,F7.3/' ','QHsteam',2X,F8.3,2X,F7.3/' ','QHtotal', 
+2X,F8.3,1X,F8.3/’ ’,'Qhrsg',4X,F8.3,2X,F7.3/’ ','Bexhaust', 
+1X,F8.3/' ','We(1)',4X,F8.3/' ',’We(2)',4X,F8.3/' ',’We(3)',4X, 
+F8.3/'0') 

IF(CYCLE.EQ.'S') THEN 

WRITE(15,185) Wtl,Bia,Biw,Be1,Wes(1), 'es(2),Wes(3),Bs1,Bs2 
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For data output of cycle efficiency and specific net work 
versus steam/air mass ratio (0.01-0.50), (3TIG cycle). 

IF(OUTPUT.EQ.'2') THEN 
Mr=0.01 

DO 270 JJ=1,491 

CALL CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,Wt,Wp,Wnet,Cph, 
+Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wt1,Bia,Biw,Be1, 
+Wes,K1 ,T5s1 ,Tstak1 ,Qhexa,T9a,Bs1 ,Bs2) 

IF(LAW.EQ.I) WRITE(17,255) Mr,ETA,Wnet,K 
FORMAT(' ',F5.5,2X,F6.4,2X,F8.3,2X,I1) 

IF(LAW.EQ.2) WRITE(17,260) Mr,(EP3(I), I=1,6),K,K1 
FORMAT(' ',F5*3,2X,6(F6.4,2X),I1,2X,I1 ) 

Mr=Mr+0.001 

CONTINUE 

ENDIF 




•V- ».’fi 
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For output of maximum First Law efficiency and corresponding mass 
ratio & specific net work versus turbine inlet temperature (STIC 
cycle). 

IF(OUTPUT.EQ.'3') THEN 
T6=1000.0 
IF(L.EQ.I) R=823 
IF(L.EQ.2) R=1050 
IF(L.EQ.3) R=1178 
IF(L.EQ.4) R=1266 
IF(L.EQ.5) R=1334 
DO 310 11=1 ,R 

IF(T6.LT.1073••AND.L.EQ.4) THEN 
Mr=0.010 
C=491 
GO TO 275 
ENDIF 

-IF(T6.LT. 1160..AND.L.EG.5) THEN 
IF(T6.LT.1066.) THEN 
Mr=0.414 

CALL CYCLE2(T6,L,Mr,CYOLE,T,P,PL,ETA,EPS,Wc,Wt,Wp. 
+Wnet,Cph,Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wtl ,Bia,Biw, 
+Be1,Wes,K1,T5s1,Tstak1 ,Qhexa,T9a,Bs1 ,Bs2) 

MAX(11,1 )=T6 
MAX(II,2)=Mr 
MAX(II,3)=ETA 
MAX(II,4)=Wnet 
MAX(11,5)=Tstack 
GO TO 295 
ENDIF 
Mr=0.010 
C=491 
GO TO 275 
ENDIF 

IP(L.EQ.I) THEN 

Mr=SL41 *(T6-1000. )-K). 13 
C=( 3142 * (T 6 - 1000.) -hj . 16-M r)/. 001 +1 . 0 


m. i 
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IF(L.EQ.2) THEN 

Mr=SL81 *(T6-1000.)+0.06 
C=(SL82*(T6-1000. )+0.10-Mr)/.001+1.0 
ENDIF 

IF(L.EQ.3) THEN 

Mr=SL121*(T6-1000.)+0.02 
C=(SL122*(T6-1000.)+0.05-Mr)/.001+1.0 
ENDIF 

IF(L.EQ.4) THEN 

Mr=SL161*(T6-1000.)+0.01 
C=(SL162*(T6-1000.)40.05-Mr)/.001+1.0 
ENDIF 

IF(L.EQ.5) THEN 

Mr=SL201*(T6-1150.)+0.025 
C=(SL202*(T6-1150.)+0.075-Mr)/.001+1.0 
ENDIF 

275 DO 280 JJ=1,C 

CALL CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,¥c,¥t,¥p,¥net, 
+Cph,Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,¥e,¥t1,Bia,Biw,Be1, 
+¥es,K1 ,T5s1 ,Tstak1,Qhexa,T9a,Bs1,Bs2) 

MAXI(JJ,1 )=Mr 
MAXI(JJ,2)=ETA 
MAXI(JJ,3)=Wnet 
MAXI (JJ,4)=Tstack 
Mr=Mr+0.001 
280 CONTINUE 

MAX(II,3)=MAX1(1,2) 

DO 290 1=2 C 

IF(MAX1(I,2).GT.MAX(II,3)) THEN 
MAX(II,1)=T6 
MAX(II,2)=MAX1(1,1) 

MAX(II,3)=MAX1(1,2) 

MAX(II,4)=MAX1(1,3) 

MAX(II,5)=MAX1(1,4) 

ENDIF 

290 CONTINUE 

295 WRITE(18,300) (MAX(II,J), J=1,5),L 

300 FORMAT(' \F6.1,2X,F5.3,2X,F6.4,2X,F8.3,2X,F9.5,2X,I1) 

T6=T6+1.0 
310 CONTINUE 
INDIF 
C 

320 STOP 
END 
C 

C-SUBROUTINES- 

C 

c 

SUBROUTINE CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,¥t,Wp,Wnet,Cph, 
+Cpec,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,¥e,¥t1,Bia,Biw,Be1,¥es, 
+K1 ,T5s1 ,Tstak1,Qhexa,T9a,Bs1 ,Bs2) 


C 

C 


This subroutine calculates the state points for the desired cycle 





















































using First Law energy balances. The work outputs/inputs and the 
heat inputs are calculated. Frcm these values, the First and 
Second Law efficiencies are calculated. 

CHARACTER* 1 CYCLE, CYCLE1 
INTEGER I,J,K,K1 ,L 

REAL T6,T(9),P(9),To5,Tn5,T5s,ETAc,ETAt,ETAp,EFA,Wc(2),Wt(2),Wp, 
+Kplost,Plost,Pr,Pro(5),PL,Mr,Cpa,Cps,Cpmix,Cph,Cpsc,Qha(2),Qhs, 
+Qh,Qhex,RH0s,Tstack,Wnet,Cpa1,Cpa2,Cpa3,Cpa4,Cpm1,Tsat1,EPS(6), 
+Bia,Biw,Be,We(3),Ya,Ys,Te1 ,Te2,Cpm2,Cpm3,Cpm4,T9a,Wt1 ,Bs1 ,Bs2,Mr1 
4Ves(3),T5s1 ,Tstak1 ,Qsat 

DATA Pro(1 )/4.0/,To5/0.0/,ETAc/0.er7/,ETAt/0.89/,ETAp/0.70/, 
+RH0s/997.0/,Kplost/12.128/,MWa/28.97/,MWs/l8.015/,Ru/8.31434/, 
+Te/0.0/ 


T(1)=300.0 

P(1 )=101 .325 

T(6)=T6 
DO 10 1=2,5 

Pro(I)=Pro(I-1)+4-0 
CONTINUE 

P(9)=P(1) 

Pr=Pro(L) 

PL=Pro(L) 

IF(CYCLE.EQ.’I') Pr=SQRT(Pr) 

Determine required compressor specific work (Wc). 

CALL CWQRK(T(1),P(1),Pr,ETAc,P(2),T(2),Wc(1)) 

IF(CYCLE.NE.'I') THEN 

T(5)=T(2) 

P(5)=P(2) 

ENDIF 

If IGT cycle is called for, find results of intercooler 
and 2d compressor. 

IF(CYCLE.EQ.’I') THEN 

Determine temp drop and pressure lost in intercooler. 

T(3)=327.5 

Plost=((T(1)+T(2))/(2.0*P(2)))*Kplost 
P(3)=P(2)-Plost 

Determine required work in 2d compressor. 

CALL CWQRK(T(3),P(3),Pr,ETAc,P(4),T(4),Wc(2)) 

ENDIF 


Set an initial guess for regenerator outlet temp, T(5)* 











































IF(CYCLE.EQ. 'I') T(5)=T(4)+100.0 

If analyzing IGT cycle, the program returns here to 
iterate the effect of the regenerator. 

Determine pressure lost across the regenerator. 

IF(CYCLE.EQ.'I') THEN 

Plost=((T(4)+T(5))/(2.0*P(4)))*Kplost 

P(5)=P(4)-Plost 

ENDIF 

Determine pressure loss across combustor. 

Plost=((T(5)+T(6))/(2.0*P(5)))*Kplost 
P(6)=P(5)-Plost 

IF(CYCLE.EQ.’I’) P(7)=P(6)/SQHT(P(6)/P(9)) 

IF(CYCLE.NE.'I') P(7)=P(9) 

Determine turbine specific work (Wt). 

CALL TWQRK(CYCLE,T(6),P(6),P(7),Mr,ETAt,T(7),Cpmix,Wt(1)) 
IF(CYCLE.NE.'I') T(9)=T(7) 

If IGT cycle is called for, find results of reheat 
and 2d turbine. 

IF(CYCLE.EQ.'I') THEN 
T(8)=T(6) 

Determine pressure loss in reheat. 

Plost=((T(7)+T(8))/(2.0*P(7)))*Kplost 
P(8)=P(7)-Plost 
Determine 2d turbine work. 

CALL TWCRK(CYCLE,T(8),P(8),P(9),Mr,ETAt,T(9),Cpmix,Vt(2)) 
ENDIF 

Determine steam exit temperature from HRSG, required heat input 
to raise exiting steam to turbine inlet temperature, and pump 
work. 

IF(CYCLE.EQ. 'S’) THEN 

CALL HRSG1 (T(1),T(6 ),T(9),Mr,L,Cph,Qhex,T5s,Cpsc,Qhs,Tstack,K 
■fTsat 1 f Qsat) 

Plost=((T(6)+T(1 ))/(2.0*P(5)))*Kplost 
Wp=Mr*(P(5)-P(1)+Plost)/(RHOs*ETAp) 

EJDIF 

Determine heat input to air in combustor. 

CALL CPAIR1(T(5)tT(6),Cpa) 















































Qha(1)=0.95*Cpa*(T(6)-T(5)) 


82 


C 

C If IGT cycle is called for, add in heat input 

C for the reheat. 

C 

IP(CYCLE.EQ.'I *) THEN 

CALL CPAIR1(T(7),T(8),Cpa) 
Qha(2)=0.95*Cpa*(T(8)-T(7)) 

ENDIP 

C 

C Iteration to account for regenerator is done here. 

C 

IF(CYCLE.EQ.'I') THEN 
Tn5=T(5) 

IP(ABS(Tn5-To5).LT.0.1) GO TO 30 
To5=Tn5 

T(5)=T(9)-50.0 
GO TO 20 
ENDIP 

C 

C Iteration to determine Texhaust from regenerator. 

C 

30 IF(CYCLE.EQ.'I') THEN 

CALL CPAIR1(T(4),T(5),Cpa1) 

Te2=T(4)+100.0 

35 CALL CPAIR1(Te2,T(9),Cpa2) 

Te2=T(9)-Cpa1*(T(5)-T(4))/Cpa2 
IF(ABS(Te2-Te1).GE.0.1) THEN 
Tel=Te2 
GO TO 35 
INDIF 
ENDIF 

C 

IF(CYCLE.NE.'I') THEN 
Wc(2)=0.0 
Wt(2)=0.0 
Qha(2)=0.0 
DO 40 J=2,4 
P(J)=0.0 
T(J)=0.0 

40 CONTINUE 

DO 50 1=7,8 
P(I)=0.0 
T(I)=0.0 

50 CONTINUE 

ENDIF 

C 

IF(CYCLE.NE.'3') THEN 
Wp=0.0 
Qhs=0.0 
Qhex=0.0 
Cph=0.0 
Cpec=0.0 
Tstack=0.0 






























ENDIF 

Determine First law cycle efficiency (ETA). 

Wnet=Wt(1)+Wt(2)-Wc(1)-¥c(2)-Wp 
ETA=Wnet/(Qha(1 )+Qha(2)+Qhs) 

Determine Second Law efficiency for desired cycle. 

IF(CYCEE.EQ.'G') THEN 

CALL CPAIR1(T(1 ),T(9),Cpa1 ) 

CALL CPAIR1(T(5),T(6),Cpa2) 

Be=0.95*Cpa1*(T(9)-T(1 )*(1 .O+L0G(T(9)/T(1)))) 

We(1)=0.95*Cpa2*(T(6)-T(5)-T(1 )*LOG(T(6)/T(5))) 

EPS(1 )=Wnet/(We(1 )-Be) 

EPS(2)=Wnet/We(1) 

ENDIF 

IF(CYCLE.EQ.'I') THEN 

CALL CPAIR1 (T(1 ),T(9),Cpa1) 

CALL CPAIR1(T(5),T(6),Cpa2) 

CALL CPAIR1 (T(7),T(8),Cpa3) 

CALL CPAIR1(T(3),T(8),Cpa4) 

Be=0.95*Cpa1 *(Te2-T(1 )*(1 .OLOG(Te2/T(1)))) 

We(1)=0.95*Cpa2*(T(6)-T(5)-T(1 )*L0G(T(6)/T(5))) 
We(2)=O.95*Cpa3*(T(8)-T(7)-T(1)*L0G(T(8)/T(7))) 
We(3)=O.95*Cpa4*(T(3)-T(2)-T(1)*L0G(T(3)/T(2))) 

EPS(1)=Wnet/(We(1)+We(2)+We(3)-Be) 

EPS(2)=Wnet/(We(1 )+We(2)+We(3)) 

ENDIF 

IF(CYCLE.EQ.'S') THEN 

CALL CFM(T(1).Tstack.Mr.Cpml ) 

CALL CPH20(T(1),T(6),Cps1) 

CALL CPH20(T(1 ),T5s,Cps2) 

CALL CPAIR1 (T(5),T(6),Cpa) 

Ysu= (0.95*MWs) / (0.95*MWs+Mr*MWa) 

Ys=(Mr*M - rfa)/(0.95*MWa-^r*MWa) 

Bia=T(1 )*Ru*LOG( 1.0/Ya)/MWa 
Biw=Mr*T (1 )*Ru*L0G(1.0/Ys)/MWs 

Be=(Mr+0.95)*Cpm1*(Tstack-T(1 )*(1 .O+L0G(Tstack/T(1)))) 

CYCLE1 ='G' 

Mr 1=0.0 

CALL TWCRK(CYCLE1 ,T(6),P(6),P(9),Mr1 ,ETAt,T9a,Cpm2,Wt1) 
We(1)=0.95*Cpa*(T(6)-T(5)-T(1)*L0G(T(6)/T(5))) 
IF(K.EQ.3.0H.K.EQ.4) THEN 
We(2) = (1.0—T(1 )/Tsat1 )*Qsat 

We(3)=Mr*Cpec*(T(6)-T5s-T(1)*L0G(T(6)/T5s)) 

CALL HRSG2(T(1),T(6),T9a,Mr,L,T5s1,Tstak1,Wes,K1,Qhexa,Bs1 ,Bs2) 
CALL CPAIR1(T(1),Tstak1,Cpa1) 

Bel=0.95*Cpa1*(Tstak1-T(1)*(1.0+LOG(Tstak1 /T(1)))) 

EPS(1)=Wnet/(We(1)4We(2)-+We(3)+Bia+Biw-Be) 
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EPS(2)=Wnet/(We(1 )+We(2)+We(3)+Bia+Biw) 

EPS(3)=(Wt1-Wc(1)-Wp+Bs1)/(We(1)+Ves(1 )+Wes(2)4¥es(3)-Be1) 
EPS(4)=(Wt1-Wc(1 )-WpfBs1 )/(We(1 )+Wes(1 )+Wes(2)+¥es(3)) 
EPS(5)=(Wt1-Wc(1 )-Wp+Bs2)/(We(1 )-Be1) 

EPS(6)=(Wt1-Wc(1)-WpfBs2)/We(1) 

ENDIF 

RETURN 

END 


SUBROUTINE CWQRK(T1,P1,Pr,ETAc,P2,T2,Wc) 

This subroutine finds the required compressor work input. 
CHARACTERS TYPE 

REAL T1,P1,Kao,Ka,Pr,P2,Ts2,ETAc,T2o,T2n,T2,Cpa,Ra,Tisen,Tact 

DATA Kao/1.397/,Ra/.28700/ 

TYPE='COMP' 

Ka=Kao 
P2=P1 *Pr 

Ts2=Tisen(TYPE,T1,Pr,Ka) 

T2o=Tact(TYPE,T1,Ts2,ETAc) 

CALL CPAIR1(T1,T2o,Cpa) 

Ka=Cpa/(Cpa-Ra) 

Ts2=Tisen(TYPE,T1 ,Pr,Ka) 

T2n=Tact(TYPE,T1,Ts2,ETAc) 

IF(ABS(T2n-T2o).GE.0.1) THEN 
T2o=T2n 
GO TO 10 
ENDIF 
T2=T2n 

Wc=Cpa*(T2-T1) 

RETURN 

END 


SUBROUTINE TW0RK(CYCLE,T1,PI ,P2,Mr,ETAt,T2,Cpmix,Wt) 

This subroutine finds the turbine work output. 

CHARACTER*! CYCLE 
CHARACTER*4 TYPE 

REAL T1,PI,P2,Mr,Kao,Ka,Kso,Ks,ETAt,Ra,Rs,T2,Cpmix,Wt,Ts2,Tisen, 
+T2o,T2n,Tact,Cpa,Cpao,Cpe,Cpso,Kmix,Pr 

DATA Kao/1-397/,Kso/l-327/,Cpao/1.0035/,Cpso/1.8723/,Ra/.28700/, 






















































+Rs/.46152/ 


TYPE='TURB' 

T2o=0.0 

Ka=Kao 

Ks=Kso 

Cpa=Cpao 

Cps=Cpso 

Pr=P1 /P2 

IF(CYCLE.NE.'S') Kmix=Ka 

IP(CYCLE.EQ.•S') Kmix=(Mr*Cps+.95*Opa)/((Mr*Cps/Ks)+(.95*Cpa/Ka)) 
Ts2=Tisen(TYPE,T1,Pr,Kmix) 

T2n=Tact(TYPE,T1,Ts2,ETAt) 

IF(ABS(T2n-T2o).GE.0.1 ) THEN 
T2o=T2n 

CALL CPAIR1(T1,T2o,Cpa) 

CALL CPH20(T1,T2o,Cps) 

Ka=Cpa/(Cpa-Ra) 

Ks=Cps/(Cps-Rs) 

GO TO 10 
END IF 
T2=T2n 

IP(CYCLE.NE.'S') Cps=0.0 
Cjmix=(Mr*Cps+0.95*Cpa)/(0.95-hMr) 

Wt=Cpmix*(.95+Mr)*(T1-T2) 

RETURN 

END 


SUBROUTINE HRSG1(T1,T6,T9,Mr,L,Cph,Qhex,T5s,Cpec,Qhs,Tstack,K, 
+Tsat1 ,Qsat) 

This subroutine calculates the performance of the HRSG for the 
STIG cycle. It finds the exiting steam temperature, the exhaust 
gas stack temperature,and the heat required to raise the exiting 
steam to the turbine inlet temperature. For this subroutine the 
hot exhaust gas is a mixture of steam and air. 

INTEGER K,L 

REAL Tsat(5),Hfg(5),Cpw(5),Q1,Q2,Q3,Qhex,Qhex1,Qhs,Tpp,Tpp1,Tpp2, 
+T1,T6,T9,Tstack,Tstak1,Tstak2,Tstak3,T5s,T5s1,T5s2,Cprnix,Cpmix1, 
+Cpmix2,Cps,Cpsc,Cph,Ch1 ,Ch2,Mr,Tsat1,Qsat 

DATA Tsat(1)/417.3/,Tsat(2)/444.1/,Tsat(3)/461.7/,Tsat(4)/475-2/, 
+Tsat(5)/486.2/,Hfg(1)/2132.4/,Hfg(2)/2046.2/,Hfg(3)/l984.0/, 
+Hfg(4)/1932.7/,Hfg(5)/l887.8/,Tpp1/0.0/,T5s1/0.0/,Tstakl/0.0/, 
-KJpw(1)/4.205/,Cpw(2)/4.215/,Cpw(3)/4.225/,Cpw(4)/4.230/, 
•+Cpw(5)/4.245/ 

Tsatl=Tsat(L) 

Check if water can reach saturation. 

























Tpp=T1+26.0 

Q1=Mr*Cpw(L)*(Tsat(L)-T1) 

CALL CPM(Tpp,T9,Mr,Cpmix) 

Cpmixl =Cpmix 

Qhexl =(Mr+0.95)*Cpmix*(T9-Tpp) 

CALL CPH20(Tsat(L),T6,Cps) 

Cpsc=Cps 

IP(Qhexl.LT.Q1) THEN 

Qhs=Q1 -Qhexl -+Mr*(Hfg(L)+Cpsc*(T6-Tsat (L))) 

Tstack=Tpp 

An approximation (due to Cpw value) of water temp entering 
combustor. 

T5s=Qhex1/(Mr*Cpw(L))+T1 
IF((T9-T5s).LT.26.0) THEN 

Pinch point is at approach temperatures. 

T5s=T9-26.0 
Tstak2=T1+50.0 
Qhex=Mr*Cpw(L)*(T5s-T1) 

IF(Qhex.GT.Q1) THEN 

T5s=Qhex1/(Mr*Cpw(L))+T1 
GO TO 15 
ENDIP 

CALL CFM(Tstak2,T9,Mr,Cpmix) 

Tstak2=T9-Qhex/((Mr+O.95)*Cpmix) 

IF(ABS(Tstak2-Tstak1).GE.0.1) THEN 
Tstakl=Tstak2 
GO TO 10 
ENDIF 

Tstack=Tstak2 

Cph=Cpmix 

Qhs=Mr*(Cpsc*(T6-Tsat(L))+Hfg(L)+Cpw(L)* 

+(Tsat(L)-T5s)) 

K=2 


The pinch point location is at the approach temperature 
location and the water exits the HRSG as a liquid. 

GO TO 90 
ENDIP 
Qhex=Qhex1 
Cph=Cpmix1 
K=1 


The pinch point location is at the water entry point and 
the water exits the HRSG as a liquid. 

GO TO 90 
ENDIP 


Check if steam can reach superheat. First assume pinch point is 
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C at entry to saturation region. 

C 

17 Tpp=Tsat(L)+26.0 

Tstakl =0.0 

20 Tstak2=Tpp-Q1 /((Mr-K).95)*Cpmix) 

IF(ABS(Tstak2-Tstak1 ).GE.0.1) THEN 
CALL CPM(Tstak2,Tpp,Mr,Cpmix) 

Tstakl=Tstak2 
GO TO 20 
ENDIF 

Tstak3=Tstak2 

CALL CPM(Tstak2,T9,Mr,Cpmix2) 

Chi =Cpmix1 *(Mr+0.95) 

Ch2=Cpmix2*(Mr+0.95) 

Qhex=Ch2*(T9-Tstak2) 

Q2=Mr*(Hfg(L)+Cpw(L)*(Tsat(L)-T1)) 

IF(Qhex.LT.Q2) THEN 
Cph=Cpmix2 

Qhs=Q2-Qhex+Mr*Cpsc*(T6-Tsat(L)) 

Tstack=Tstak2 
Qsat=Q2-Qhex 
K=3 

The pinch point location is at the beginning of the sat¬ 
uration region and the steam exits the HRSG in saturation. 

Check to see if pinch point is at water entry point. 

IF((Tstack-TI).LT.26.0) THEN 
Tstack=T1+26.0 

Check to see if Mode 3 is more applicable. 

Tpp2=Tsat(L)+30.0 
CALL CPM(Tstack,Tpp2 f Mr,Cpmix) 

Tpp2=Q1 /((Mr+0.95)'*'Cpmix)+Istack 
IF(ABS(Tpp2-Tpp1 ).GE.0.1) THEN 
Tppl =Tpp2 
GO TO 25 
ENDIF 

IF((Tpp2-Tsat(L)).LT.26.0) THEN 
GO TO 17 
ENDIF 

Qhex=Ch1*(T9-Tstack) 

Cph=Cpmix1 

Qhs=Q2-Qhex+Mr*Cpsc*(T6-Tsat(L)) 

Qsat=Q2-Qhex 
K=4 

The pinch point location is at the water entry point and the 
steam exits the HRSG in saturation. 



















MDIF 
GO TO 90 
MDIF 

last choice, steam reaches superheat with pinch point at water 
entry point, entry to saturation region, or at steam exit 
(approach temperatures). 

With assumption pinch point is at approach temperatures check to 
see if temperature difference at other possible pinch point loca¬ 
tions is more than minimum limit. 

T5s=T9-26.0 

CALL CPH20(Tsat(L),T5s,Cps) 

Q3=Mr*(Cps*(T5s-Tsat(L))+Hfg(L)) 

Initial guess on Tpp2 to determine Cpmix. 

Tppl=0.0 

Tpp2=Tsat(L)+100.0 
CALL CPM(Tpp2,T9,Mr,Cpmix) 

Tpp2=T9-Q3/((.95+Mr)*Cpmix) 

IF(ABS(Tpp2-Tpp1 ).GE.0.1) THEN 
Tppl=Tpp2 
GO TO 30 
MDIF 
Tpp=Tpp2 

Qhex=Q3+Q1 
Tstakl=0.0 

Tstak2=T9-Qhex/(Cpmix*(Mr+0.95)) 

IF(ABS(Tstak2-Tstak1 ).GE.0.1) THM 
CALL CIW(Tstak2,T9,Mr,Cpmix) 

Tstakl =Tstak2 
GO TO 40 
ENDIF 

Tstack=Tstak2 

Cph=Cpmix 

CALL CPH20(T5s,T6,Cpsc) 

Qhs=Mr*Cpec*(T6-T5s) 

K=5 

The pinch point location is at the approach point temperature 
location and the steam exits the HRSG as a superheated vapor. 

Check to see if Mode 6 or 7 are more applicable. 

IF((Tstack-T1).LT.26.0.AND.(Tstack-TI).LT.(Tpp-Tsat(L))) GO TO 70 
IF((Tpp-Tsat(L)).LT.26.0) THM 
Tpp=Tsat(L)+26.0 
Tstack=Tstak3 
Cph=Cpmix2 
Qhex=Ch2*(T9-Tstack) 

T5s2=Tsat(L)+100.0 
























CALL CPH20(Tsat(L),T5s2,Cps) 

T5s 2= (Qhex-Q2) / (Mr*Cps )+Tsat (L) 

IF(ABS(T5s2-T5s1).GE.0.1) THEN 
T5s1=T5s2 
GO TO 60 
ENDIF 
T5s=T5s2 

CALL CFH20(T5s,T6,Cpsc) 

Qhs=Mr*Cpsc*(T6-T5sj 

IF((Tstack-TI).LT.26.0) THEN 
GO TO 70 
ENDIF 

K=6 

The pinch point location is at the beginning of the sat¬ 
uration region and the steam exits as a superheated vapor 

GO TO 90 
ENDIF 

GO TO 90 

Tstack=T1+26.0 

Check to see if Mode 6 is more applicable. 

Tpp1=0.0 

Tpp2=Tsat(L)+30.0 

CALL CFM(Tstack f Tpp2,Mr,Cpmix) 

Tpp2=Q1 /((Mr+0.95)*Cpmix)+Tstaok 
IF(ABS(Tpp2-Tpp1 J.GE.0.1) THEN 
Tppl=Tpp2 
GO TO 75 
ENDIF 

IF((Tpp2-Tsat(L)).LT.26.0) THEN 
GO TO 50 
ENDIF 

CALL CPM(Tstack,T9,Mr,Cph) 

Qhex=(Mr+O.95)*Cph*(T9-Tstack) 

T5s1=0.0 

T5s2=T9-50.0 

CALL CFH20(Tsat(L),T5s2,Cps) 

T5s2=(Qhex-Q2)/(Mr*Cps)+Tsat(L) 

IF(ABS(T5s2-T5s1J.GE.0.1) THEN 
T5s1=T5s2 
GO TO 80 
ENDIF 
T5s=T5s2 

CALL CFH20(T5s,T6,Cpsc) 

Qhs=Mr*Cpec*(T6-T5s) 









































The pinch point location is at the water entry point and the 
steam exits the HRSG as a superheated vapor. 


RETURN 

END 


SUBROUTINE HRSG2(T1,T6,T9,Mr,L,T5s,Tstack,Wes,K1 ,Qhex,Bs1 ,Bs2) 

This subroutine calculates the performance of a HRSG combined 
with a GT cycle to produce process steam. In this case the hot 
exhaust gas is air only as there is no steam injection. This 
subroutine does the same calculations as the HRSG1 subroutine. 

INTEGER K1 ,L 

REAL T1,T6,T9,Tpp,Tpp1,Tpp2,Tstack,Tstak1,Tstak2,Tstak3,T5s, 
+T5s1 ,T5s2,Wes(3),Cpa,Cpa1 ,Cpa2,Cps,Cps1,Ch1,Ch2,Q1 ,Q2,Q3,Qhex, 
+Qhex1 ,Tsat(5),Hfg(5),Cpw(5),Sfg(5),Mr,Bs1,Bs2,X 

DATA Tstakl/0.0/,Tppl/0.0/,T5sl/0.0/,Tsat(1)/417-3/, 
+Tsat(2)/444.1/,Tsat(3)/46l.7/,Tsat(4)/475-2/,Tsat(5)/486.2/, 
4flfg(l)/2132.4/,Hfg(2)/2046.2/,Hfg(3)/l984.0/,Hfg(4)/1932.7/, 
4Hfg(5)/l887.8/,Cpw(l)/4.205/,CpwT2)/4.215/,Cpw(3)/4.225/, 
-fOpw(4)/4.230/,Cpw(5)/4.245/,Sfg(1 )/5 • 1100/, Sfg(2)/4.6067/, 

4Sfg(3)/4.2964/,Sfg(4)/4•0670/,Sfg(5)/3•8825/ 

CALL CPH20(Tsat(L),T6,Cps) 

Bsl =Mr*(Cpw(L)*(Tsat(L)-T1 *(1 .O+L0G(Tsat(L)/T1 )))+ 

+Hfg(L)-T1 *Sfg(L)+Cps*(T6-Tsat (L)-T1 *L0G(T6/Tsat (L)))) 

Tpp=T1+26.0 

Q1 =Mr*Cpw(L)*(Tsat(L)-T1) 

CALL CPAIR1 (Tpp,T9,Cpa) 

Cpal=Cpa 

Qhexl =0.95*Cpa1*(T9-Tpp) 

CALL CPH20(Tsat(L),T6,Cps1 ) 

IF(Qhexl.LT.Q1) THEN 
Tstack=Tpp 

T5s=Qhex1/(Mr*Cpw(L))+T1 

Wes(1)=Mr*Cpw(L)*(Tsat(L)-T5s-T1*L0G(Tsat(L)/T5s)) 
Wes(2)=Mr*(1.0-T1 /Tsat(L))*Hfg(L) 

Wes(3)=Mr*Cps1*(T6-Tsat(L)-T1*L0G(T6/Tsat(L))) 
Bs2=Mr*Cpw(L)*(T5s-T1*(1 .Q+L0G(T5s/T1 ))) 

K1 =1 

IF((T9-T5s).LT.26.0) THEN 
T5s=T9-26.0 

Qhex=Mr*Cpw(L)*(T5s-T1) 

IF(Qhex.GT.Q1 ) THEN 

T5s=Qhex1/(Mr*Cpw(L))+T1 
GO TO 20 
ENDIF 

Tstak2=T1+50.0 

CALL CPAIR1(Tstak2,T9,Cpa) 

Tstak2=T9-Qhex/(0.95*Cpa) 

























































IF(ABS(Tstak2-Tstak1).GE.0.1) THEN 
Tstakl=Tstak2 
GO TO 10 
ENDIF 

Ts“tack=Ts'tak2 

Wes(1)=Mr*Cpw(L)*(Tsat(L)-T5s-T1*LOG(Tsat(L)/T5s)) 
Wes(2)=Mr*(1.0-T1 /Tsat(L))*Hfg(L) 

Wes(3)=Mr*Cps1*(T6-Tsat(L)-T1*L0G(T6/Tsat(L))) 
Bs2=Mr*Cpw(L)*(T5s-T1 *(1.0+LOG(T5s/T1))) 

K1 =2 

GO TO 130 
ENDIF 
Qhex=Qhex1 
GO TO 130 
ENDIF 

Tpp=Tsat(L)+26.0 

fPo + olfl —O O 

Tstak2=Tpp-Q1 /(0.95*Cpa) 

IF(ABS(Tstak2-Tstak1 ).GE.0.1) THEN 
CALL CPAIR1(Tstak2,Tpp,Cpa) 

Tstakl =Tstak2 
GO TO 40 
ENDIF 

Tstak3=Tstak2 

CALL CPAIR1 (Tstak2,T9,Cpa2) 

Chi =Cpa1 *0.95 
Ch2=Cpa2*0.95 
Qhex=Ch2*(T9-Tstak2) 

Q2=Mr*(Hfg(L)+Cpw(L)*(Tsat(L)-T1)) 

IF(Qhex.LT.Q2.QR.Qhex1 .LT.Q2) THEN 
T5s=Tsat(L) 

Wes(2)=(1•0-T1/Teat(L))*(Q2-Qhex) 

Wes(3)=Mr*Cps1 *(T6-T5s-T1 *L0G(T6/T5s)) 

X=(Qhex-QI)/(Mr*Hfg(L)) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1.0+LOG(Tsat(L)/T1)))- 
+T1 *X*Sfg(L))+Qhex-Q1 
K1 =3 

IF((Tstack-T1).LT.26.0) THEN 
Tstack=T1+26.0 
Tpp2=Tsat(L)+30.0 
CALL CPAIR1(Tstack,Tpp2,Cpa) 

Tpp2=Q1/(0.95*Cpa)+Tstack 
IF(ABS(Tpp2-Tpp1 ).GE.0.1) THEN 
Tppl=Tpp2 
GO TO 50 
ENDIF 

IF((Tpp2-Tsat(L)).LF.26.0) GO TO 30 
Qhex=Ch1*(T9-Tstack) 

Wes(2)=(1.0-T1/Tsat(L))*(Q2-Qhex) 

Wes(3)=Mr*Cps1*(T6-T5s-T1 *LOG(T6/T5s)) 

X=(Qhex-QI )/(Mr*Hfg(L)) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1.0+L0G(Tsat(L)/T1 )))- 
+T1 *X*Sfg(L))+Qhex-Q1 










































GO TO 130 
MDIF 

T5s=T9-26.0 

CALL CPH20(Tsat(L),T5s,Cps) 
Q3=Mr*(Cps*(T5s-Tsat(L))+Hfg(L)) 

Tpp1=0.0 

Tpp2=Tsat(L)+100.0 
CALL CPAIR1(Tpp2,T9»Cpa) 

Tpp2=T9-Q3/(0.95*Cpa) 

IF(ABS(Tpp2-Tpp1J.SE.0.1) THEN 
Tppl=Tpp2 
GO TO 60 
ENDIF 
Tpp=Tpp2 
Qhex=Q3+Q1 
Tstak1=0.0 

Tstak2=T9-Qhex/(Cpa*0.95) 

IF(ABS(Tstak2-Tstak1j.GE.0.1) THEN 
CALL CPAIR1(Tstak2,T9,Cpa) 

Tstakl=Tstak2 
GO TO 70 
INDIP 

Tstack=Tstak2 

CALL CFH20(T5s,T6,Cps) 

Wes(3)=Mr*Cps*(T6-T5s-T1*L0G(T6/T5s)) 

CALL CPH20(T5s,Tsat(L),Cps) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1 .OfLOG(Tsat(L)/T1)))+ 
•+Hfg(L)-T1 *Sfg(L)40ps*(T5s-Tsat(L)-T1*L0G(T5s/Tsat(L)))) 
K1 =5 

IF((Tstack-TI).LT.26.0.AND.(Tstack-TI).LT.(Tpp-Tsat(L))) 
-KX) TO 100 

IF((Tpp-Tsat(L)).LT.26.0) THEN 
Tpp=Tsat(L)+26.0 
Tstack=Tstak3 
Qhex=Ch2*(T9-Tstack) 

T5s2=Tsat(L)+100.0 

CALL CPH20(Tsat(L),T5s2,Cps) 

T5s2=(Qhex-Q2)/(Mr*Cps)+Tsat(l) 

IF(ABS(T5s2-T5s1 J.GE.0.1) THEN 
T5s1=T5s2 
GO TO 90 
MDIF 
T5s=T5s2 

CALL CPH20(T5s,T6,Cps) 

¥ea(3)=Mr*Cps*(T6-T5s-T1 *L0G(T6/T5s)) 

CALL CPH20(Tsat{L),T5s.Cpe) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1 *(1 .OfLOG(Tsat(L)/T1 )))+ 
4Hfg(L)-T1*Sfg(L)4Cpe*(T5s-Tsat(L)-T1*L0G(T5s/Tsat(L)))) 
K1=6 

IF((Tstack-TI).LT.26.0) GO TO 100 
GO TO 130 


























































ooooooo ooo 


93 


GO TO 130 

100 Tstack=T 1+26.0 

Tpp1=0.0 

Tpp2=Tsat(L)+30.0 
110 CALL CPAIR1 (Tstack,Tpp2,Cpa) 

Tpp2=Q1/(0.95*Cpa)+Tstack 
IF(ABS(Tpp2-Tpp1J.GE.0.1) THEN 
Tppl=Tpp2 
GO TO 110 
ENDIF 

IF((Tpp2-Tsat(L)).LT.26.0) GO TO 80 
Qhex=Ch1 *(T9-Tstack) 

T5s1=0.0 

T5s2=T9-50.0 

120 CALL CPH20(Tsat(L),T5s2,Cps) 

T5s2=(Qhex-Q2)/(Mr*Cps)+Tsat(L) 

IF(ABS(T5s2-T5s1).GE.0.1) THEN 
T5s1=T5s2 
GO TO 120 
ENDIF 
T5s=T5s2 

CALL CPH20(T5s,T6,Cps) 

Wes(3)=Mr*Cps*(T6-T5s-T1*L0G(T6/T5s)) 

CALL CPH20(Tsat(L),T5s,Cps) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1.O+L0G(Tsat(L)/T1)))+ 
+Hfg(L)-T1 *Sfg(L)+Ops*(T5s-Tsat(I)-T1 *L0G(T5s/Tsat (L)))) 
K1 =7 

130 RETURN 
END 


SUBROUTINE CPAIR1(T1,T2,Cpa) 

This subroutine finds the specific heat at constant pressure 
of air as a function of average temperature of the process. 

It is based on a relation by Irvine and Liley, STEAM AND GAS 
TABLES WITH COMPUTER EQUATIONS, Academic Press, Inc., Orlando, 
EL, 1984. 

INTEGER I 

REAL T1,T2,Tavg,A(5),Cpa 
C 

DATA A(1)/0.103409E1/,A(2)/-0.2848870E-3/,A(3)/0.7816818E-6/, 
+A(4)/-0•4970786E-9/,A(5)/0.1077024E-12/ 

C 

Cpa=0.0 

Tavg=(T1+T2)/2.0 
DO 10 1=1,5 

Cpa=Cpa+A(I)*(Tavg**(1-1)) 

10 CONTINUE 
RETURN 
END 
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SUBROUTINE CHi20(T1,T2,Cps) 

This subroutine finds the specific heat at constant pressure 
of steam (ideal gas) as a function of average temperature of 
the process. It is based on a relation by Reynolds, THERMO¬ 
DYNAMIC PROPERTIES IN SI, Stanford University, Stanford, CA, 
1979- 

INTEGER I 

REAL T1,T2,Tavg,G(6),Rs,Cvs,Cps 

DATA G(1 )/4.6E+04/,G(2)/l .011249*>+03/,G(3)/8.3893E-01 /, 

-KJ (4)/-2.19989B-04/,G(5)/2.46619E-07/,G(6)/-9.7047E-11/, 
+Rs/.46152/ 

Cvs=0.0 

Tavg=(T1+T2)/2.0 
DO 10 1=1,6 

Cvs=Cvs+G(I)*(Tavg**(1-2)) 

CONTINUE 

Cps=(Cvs/l000.)+Rs 

RETURN 

END 


SUBROUTINE CiW(T1 ,T2,Mr,Cpmix) 

This subroutine finds the specific heat at constant pressure 
of the steam/air mixture of the STIG cycle using the two 
subroutines listed above. 

REAL T1,T2,Mr,Cpa,Cps,Cpmix 

CALL CPAIR1(T1,T2,Cpa) 

CALL CPH20(T1,T2,Cps) 

Cpmix= (Mr*Cps-K). 95*Cpa) / (0.95+Mr) 

RETURN 

END 

-REAL FUNCTIONS- 

REAL FUNCTION Tisen(TYPE,T1,Pr,Ka) 

This function finds the end state temperature for an isentropic 
compression or expansion process. 

CHARACTERS TYPE 
REAL T1,Pr,Ka 



IF(TYPE.EQ.'COMP’) THEN 






























































Tisen=T1*(Pr**((Ka-1.0)/Ka)) 
TiTSP 

Tisen=T1/(Pr**((Ka-1.0)/Ka)) 
ENDIP 
RETURN 
END 


REAL FUNCTION Tact(TYPE,T1,Ts2,ETA) 

This function finds the actual end state temperature for a 
compression or expansion process using the isentropic 
process end temperature and the component's isentropic 
efficiency. 

CHARACTERS TYPE 
REAL T1,Ts2,ETA 

IF(TYPE.EQ.'COMP') THEN 
Tact=T1+(Ts2-T1)/ETA 
USE 

Tact=T1-ETA*(T1-Ts2) 

ENDIP 

RETURN 
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PROGRAM CONDEN 

This program modifies data files to insure no steam condenses 
inside the components of the STIG cycle. It checks the output 
files of max First Law efficiency and corresponding steam/air 
mass ratio versus turbine inlet temperature. If the stack temp 
at max efficiency is below the dewpoint temp then the mass ratio 
is decreased until the stack temp is above dewpoint. This program 
must be linked with the subroutines of PROGRAM GT. 

CHARACTER*! CYCLE 
CHARACTERS 3 FNAME 
INTEGER I,J,K,K1 ,L,NN,SIZE 

REAL MAX(1334,5),T(9),T6,Tstack,Tstak1,T5s,T5s1,Tdp,T9a,P(9),PL, 
+ETA,EPS(6),Wc(2),Wt(2),Wp,Wnet,Cph,Cpsc,Cpmix,Qha(2),Qhs,Qhex, 
+Qhexa,Be,We(3),Wt1,Bia,Biw,Be1,Wes(3),Bs1,Bs2,Mr 

CYCLERS' 

WRITE(*,10) 

FORMAT('0', 'What is the file to be reviewed (FO R .DAT; )?’) 
READ(*,20) FNAME 
FORMAT(A13) 

OPPN (UNIT=10,NAME=FHAME,STATUS^'UNKNOWN•) 


NN=1 

READ (10,*,END=40) (MAX(NN,J), J=1,5),L 

NN=NN+1 

GO TO 30 

CONTINUE 

CLOSE(UNIT=10,STATUS='SAVE') 

SIZE=NN-1 
DO 80 1=1,SIZE 
T6=MAX(I,1) 

Mr=MAX(I,2) 

Tstack=MAX(1,5) 

CALL DEWPT(Mr,Tdp) 

IF(Tstack.LI.Tdp) THEN 


DO 50 J=1,500 
Mr=Mr-0.001 
CALL DEWPT(Mr,Tdp) 

CALL CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,Wt,Wp, 
+Wnet,Cph,Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wt1 ,Bia, 
+Biw,Be1,Wes,K1,T5s1,Tstak1,Qhexa,T9a,Bs1 ,Bs2) 
IF(Tstack.GT.Tdp) GO TO 60 
CONTINUE 

MAX(I,2)=Mr 

MAX(I,3)=ETA 


















MAX(I,4)=Wnet 

MAX(I,5)=Tstack 

ENDIF 

WRITE(19,70) (MAX(I,J), J=1,5) 

FORMAT(' , ,F6.1,2X,F7.3,2X,P6.4,kX,F8.3,2X,F9.5) 
CONTINUE 

STOP 

END 


■SUBROUTINE 


SUBROUTINE DEWPT(Mr,Tdp) 

This subroutine determines steam saturation temperature as a 
function of saturation pressure. It is based on a relation 
by Irvine and Liley, STEAM AND GAS TABLES WITH COMPUTER 
EQUATIONS, Academic Press, Inc., Orlando, PL, 1984. 

REAL Mr,Psat,Ys,Tdp 

DATA MWa/28.97/,MWs/18.1015/,A/42.6776/,B/-3892.70/,C/-9.48654/ 

Ys=(Mr*MWa)/(0.95*MWs+Mr*MWa) 

Psat=Ys*101.325 

Tdp=A+B/(ALOG(Psat/1000.0)4C) 


RETURN 
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